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Chanter  1:  Introduction  and  Project  Objective 


Project  Objective 

The  objective  of  this  project  was  to  develop  a  model  that  can  be  used  to 
evaluate  the  desi^  and  performance  of  a  multiple  staged  anoxic-aCTobic 
wastewater  treatment  system.  This  type  of  i^stem  has  general  qrplication,  and 
was  first  used  at  the  91st  Ave.  wastewater  treatment  plant  in  Phoenix,  Arizona 
The  model  simulates  the  concentrations  of  ammonia  nitrate,  and  nitrite  in  each  of 
the  stages  of  the  ten-stage  system,  using  a  Nficrosoft  Excel  spreadsheet.  The 
oxygen  consumption  rate  is  also  computed  in  each  stage.  The  staged  treatment 
system  performance  was  evaluated  as  a  fimction  of  various  design  and  operating 
parameters  which  include  the  effect  of  Solids  Retention  Time  (SRT),  the  internal 
recycle  flow  rate  and  distribution,  and  aeration  tank  DO  concentrations.  Two 
models  were  developed  for  comparison.  The  first  model  was  based  on  hfichaelis- 
Menton  kinetics.  The  second  model  was  a  modified  version  of  the  first  model  that 
takes  into  account  soluble  substrate  uptake  and  storage  in  die  initial  stages. 


Background  Informatjon 

Nitrification  and  Denitrification  Processes 

Anoxic/aerobic  treatment  systems  are  used  in  wastewater  treatment  for  die 
control  and  removal  of  nitrogen.  Nitrogen  concentrations  in  treatment  plant 
efQuents  must  be  controlled  in  order  to  avoid  adverse  effects  in  receiving  waters. 
High  effluent  ammonia-nitrogen  concentrations  may  result  in  depletion  of 
dissolved  oxygen  in  receiving  waters,  thereby  having  a  negative  impact  on  aquatic 
life.  Significant  nitrogen  concentrations  in  effluent  can  also  accelerate  the 
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eutrophication  of  lakes  and  allow  for  the  growth  of  rooted  aquatic  plants  and  algae 
in  streams.  Therefore,  nitrogen  control  strategies  are  an  important  part  of 
wastewater  treatment  (Metcalf  and  Eddy,  Inc.,  1991 ) 

The  biological  removal  of  nitrogen  from  wastewater  requires  a  two  step 
process.  The  influent  ammonia  is  first  converted  to  nitrate  through  ^  process 
known  as  biological  nitrification.  The  transformation  of  ammonia-nitrogen  to 
nitrate-nitrogen  does  not  result  in  nitrogen  removal,  but  it  does  eliminate  its 
oxygen  demand  (Metcalf  and  Eddy,  fric.,  1991).  The  nitrate  can  then  be  used  as 
an  electron  acceptor  under  anoxic  conditions  (without  oxygen),  and  can 
subsequently  be  reduced  to  nitrogen  gas  through  a  process  known  as 
denitrification  (Stensel,  1992(a)).  When  describing  the  denitrification  process,  the 
term  "anoxic"  is  used  rather  than  "anaerobic".  This  is  because  the  principal 
metabolic  pathways  involved  in  denitrification  are  modified  aerobic  pathways 
using  nitrate  as  the  final  electron  acceptor,  rather  than  anaerobic  reactions. 
(Metcalf  and  Eddy,  Inc.,  1991). 

Nitrification  is  a  two  step  process  involving  two  different  genera  of  nitrifying 
bacteria.  In  the  first  step,  Nitrosomonas  converts  ammonia  to  nitrite.  In  the 
second  step,  nitrite  is  converted  to  nitrate  by  Nitrobacter  (Metcalf  and  Edcfy,  Inc., 
1991).  This  two  step  oxidation  process  is  generally  written  as  follows: 

2NH4'*'  +  302  ===>  2N02"  +  4H'^  +  2H20  via  Nitrosomonas 

2N02’  +  O2  ====>  2N03"  via  Nitrobacter 

Nitrification  is  an  autotrophic  process.  Ammonia  oxidation  provides  the  energy 
for  growth  of  the  nitrifying  microorganism  and  carbon  dioxide  is  used  for  cell 
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synthesis.  This  can  be  contrasted  with  heterotrophic  organisms,  which  use  organic 
carbon  for  cell  synthesis  (Metcalf  and  Eddy,  Inc.,  1991).  Dissolved  oxygen 
concentrations  above  1  mg/1  are  required  to  sustain  reasonable  nitrification  rates. 

The  anoxic  process  of  denitrification  can  be  accomplished  by  several  different 
genera  of  heterotrophic  bacteria,  such  as  Pseudomonas,  Bacillus,  Achromobacter, 
Flavobacterium,  and  Spirillum.  These  heterotrophs  are  capable  of  dissimilatory 
nitrate  reduction  in  the  form  of  the  following  reaction  sequence: 

NO3'  =>N02-  =>N0  =>.N20  =>  N2 

This  reaction  sequence  shows  that  the  first  step  of  denitrification  involves  the 
conversion  of  nitrate  to  nitrite.  This  is  followed  by  the  production  of  nitric  oxide, 
nitrous  oxide,  and  nitrogen  gas.  These  three  gaseous  products  can  be  released  to 
the  atmosphere,  thereby  completing  the  nitrogen  removal  process  (Metcalf  and 
Eddy,  Inc.,  1991).  The  presence  of  dissolved  oxygen  wiU  suppress  the  enzyme 
system  needed  for  this  dissimilatory  nitrate  reduction  process  to  occur. 


Phoenix  91st  Avenue  Anoxic-Aerobic  Wastewater  Treatment  Plant 


The  City  of  Phoenix  91st  Avenue  Wastewater  Treatment  Plant  provides 
service  for  seven  cities  in  the  Phoenix,  AZ  area,  receiving  wastewater  that  is 
primarily  of  domestic  origin.  This  plant  uses  conventional  primary  and 
secondary  treatment  with  a  design  capacity  of  approximately  1 50  MOD.  In 
September  1992,  two  existing  aeration  tanks  were  modified  to  accomplish  nitrogen 
removal  by  converging  to  an  anoxic-aerobic  system  with  internal  recycle.  The 
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purpose  of  the  modification  was  to  achieve  biological  nitrification  and 
denitrification  without  tire  added  expense  of  constructing  new  aeration  basins. 
The  nitrogen  removal  goal  was  an  effluent  total  nitrogen  concentration  of  10.0 
mg/1  or  less.  Test  results  over  a  six  month  period  indicated  that  this  goal  was 
being  achieved  (Stensel,  1993(b)). 


Figure  1  shows  the  layout  of  one  of  the  two  modified  aeration  tanks.  The  tank 
volume  is  3.72  million  gallons  (315ftxlOOftxl5.5ft  deep),  and  it  is  sub-divided 
into  4-pass  aeration  channels  with  a  25  ft.  width.  The  total  flow  length  of  the  4- 
pass  channel  system  is  1260  ft.  The  primary  effluent  feed  rate  into  stage  1  (Ax-1) 
of  the  basin  is  approximately  20  MGD.  The  internal  recycle  flow  rate  into  stage  1 
(Ax-1)  and  stage  3  (Ax-3)  is  20  MGD  and  40  MGD  respectively.  The  return 
recycle  sludge  flow  rate  is  ^rproximately  9  MGD,  and  the  operating  solids 


retention  time  (SRT)  is  around  5  days  (Stensel,  1993(b)). 


The  basin  was  sub-divided  into  several  small  initial  anoxic  stages,  instead  of 
using  only  one  single  anoxic  stage  of  greater  volume.  This  design  feature  was 
included  in  an  effort  to  reduce  sludge  volume  index  (SVI)  values,  and  to  achieve 
shorter  total  anoxic  detention  times.  These  initial  anoxic  stages  (Ax-1  and  Ax-2), 
with  short  detention  times  and  high  F/M  ratios,  serve  as  biological  selector  zones 
that  maximize  soluble  substrate  uptake  by  non-filamentous  bacteria  (Stensel, 

1 993(b)).  This  helps  to  control  sludge  settling  characteristics.  The  first  two 
anoxic  stages  in  the  basin  (Ax-1  and  Ax-2)  are  mixed  with  course  bubble  diffusers 
The  remainder  of  the  anoxic  and  aerobic  stages  are  mixed  with  fine  bubble 
diffusers.  The  operating  dissolved  oxygen  concentration  in  the  first  four  anoxic 


zones  is  maintained  at  or  near  zero.  The  fifth  stage  (Ax-5)  is  a  very  short  zone 


with  limited  DO  concentration.  The  last  five  aerobic  zones  (Ox-2  through  Ox-6) 
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are  operated  with  target  DO  concentrations  of  2  to  3  mg/1.  SignilBlcant  biological 
storage  of  SCOD  occurs  in  the  first  stage  of  treatment  (Stensel,  1993(b)). 
Approximately  70  to  80  percent  of  the  SCOD  removed  in  the  first  anoxic  zone  is 
taken  into  cell  storage  by  the  microorganism  population  in  die  wastevs^ter. 
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FIGURE  1:  Schematic  of  Phoenix  91st  Ave. 

Anoxic-Aerobic  Treatment  System 
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Chapter  2;  Model  Based  on  Michaeiis-Menton  Kinetics 

A  Microsoft  Excel  (version  4.0)  spreadsheet  has  been  used  to  solve  the 
anoxic-aerobic  wastewater  treatment  system  model.  Ammonia,  nitrate,  and  nitrite 
concentrations  are  calculated  and  plotted  through  each  stage  of  the  ten-stage 
system.  Other  model  parameters  include  soluble  chemical  oxygen  demand 
(SCOD),  particulate  chemical  oxygen  demand  (PCOD),  and  the  oxygen 
consumption  rate  in  each  stage.  Mass  balance  equations  were  written  for 
ammonia,  nitrate,  nitrite,  SCOD  and  PCOD  for  each  of  the  ten  aeration  basin 
stages.  The  oxygen  consumption  was  based  on  the  amount  of  nitrogen  and  COD 
oxidized  in  each  stage,  and  ^ropiiate  stoichiometry  relationships  between 
oxygen  consumed  and  unit  substrate  removed  Steady-state  conditions  were 
assumed  for  the  model  solution. 

Internal  recycle  flow  is  withdrawn  fi-om  stage  9  of  ftie  aeration  basin. 
Therefore,  the  total  flow  into  the  first  anoxic  stage  includes: 

(Qir(i))  Internal  recycle  flow  fi-om  stage  9 

(Qras)  Return  flow  from  stage  10  (via  secondary  clarifier) 

(Qo)  Influent  flow  (fi-om  primary  effluent) 

The  total  flow  into  stage  3  includes  the  above  listed  flows  and  the  internal  recycle 
flow  jfrom  stage  9  (Qir(3)). 

The  following  is  a  listing  of  typical  mass  balance  equations  in  stages  1  and  6 
for  each  of  me  wastewater  components  being  modeled  Descriptions  of  kinetic 
terms  within  the  equations  are  also  included  (see  Table  1  for  coefficient  and 
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Symbol 

Definition 

Typical 

Cell# 

Value 

on 

Spread¬ 

sheet 

Mm(NS) 


Mm(NB) 


Mm(H) 


kN(ns) 


KN(nb) 


Ss(10/Xh 


maximum  specific  growth  rate  for 
nitrosomonas  nitrifiers 


maximum  specific  growth  rate  for 
nitrobacter  nitrifiers 


maximum  specific  growth  rate  for 
heterotrophs 


yield  of  nitrosomonas  nitrifiers 


yield  of  nitrobacter  nitrifiers 


yield  of  heterotrophs 


heterotrohic  maximum  specific  substrate 
utilization  rate 


half-saturation  constant  for  soluble 
substrate  (for  heterotrophs) 


DO  half-saturation  constant  for 
heterotrophs 


coefficient  for  particulate  substrate 
utilization 


DO  nitrate  reduction  inhibition  constant 


DO  half-saturation  constant  for  nitrifiers 


Nitrogen  half-saturation  constant  for 
nitrosomonas 


Nitrogen  half-saturation  constant  for 
nitrobacter 


coefficient  for  stored  substrate  utilization 


half-saturation  constant  for  stored 
substrate  utilization 


ratio  of  SCOD  in  cell  storage  to 
heterotrophic  biomass  (stage  10) 


0.50  d” 


8gCOD/g  N24 
TSS-d 


20  mg/1  N13 


0.1  mg/1 


Ss(9yXH 


Fdn 


NOo 


NO 


NO 


N02, 


N02 


N02 


ratio  of  SCOD  in  cell  stor^e  to 
heterotrophic  biomass  (stage  9) 


endogenous  decay  coef  for  heterotrophs 


endogenous  decay  coefficient  for 
nitrifiers 


nitrate-nitrogen  half-saturation  constant 


fraction  of  heterotrophs  using  nitrate  for 
electron  acceptor 


fraction  of 


dissolved  oxygen  concentration 


effluent  soluble  substrate  concentration 
for  stage  10 


soluble  substrate  concentration  in 


residual  (non-biodegradable)  SCOD 


influent  nitrate  concentration  (primary 
effluent) 


nitrate  nitrogen  concentration  in 


nitrate  concentration  in 


influent  nitrite  concentration  (primary 
effluent) 


nitrite  concentration  in  s 


nitrite  concentration  in 


influent  particulate  COD  concentration 
(primary  effluent) 


particulate  COD  concentration  in  stage 
10 


articulate  COD  concentration  in 


ammonia  concentration  in 


ammonia  concentration  in  s 


volume  of  stage  1 


volume  of  stage  2 


volume  of  stage  3 


volume  of  stage  4 


volume  of  stage  5 


volume  of  stage  6 


volume  of  stage  7 


0.06 


0.01  g/g-d 


0.1 


0/3  mg/1 


228  mg/1 


0.1  m 


30 


Omg/1 


T5 


NlSand 

N26 


N3 


6.3 


0.156  MG 


0.145  MG 


0.5047  MG  N28 


0.124  MG  |N29 


0.124  MG 


0.3394  MG 


0.4648  MG 


Vg 

volume  of  stage  8 

0.9296  MG 

N33 

V9 _ 

volume  of  stage  9 

0.4648  MG 

N34 

volume  of  stage  10 

0.4648  MG 

N35 

bBS 

total  volume  of  anoxic-aerobic  basin 

3.72  MG 

P14 

0. 

SRT  (solids  retention  time) 

5  days 

L2 

Qo 

influent  flow  (primary  effluent) 

20MGD 

N2 

Qras _ 

9.35  MGD 

N4 

internal  recycle  flow  to  stage  1 

20MGD 

N6 

internal  recycle  flow  to  stage  3 

40  MGD 

N25 

xhI _ 

heterotroph  concentration 

N12 

nitrobacter  concentration 

P4 

nitrosomonas  concentration 

T4 

Xio 

influent  inert  solids  concentration 

R6 

Xi 

inert  solids  concentration 

R7 

XT 

total  MLSS 

WJUmIiM 

P21 

Xe _ 

particulate  solids  concentration 

N/A 

** 


data  (Stensd,  H.D.,  1993(c)).  The  value  for  was  derived  by  using  a  ratio  of 
(Mm(Kb)f  t^rtamed  from  an  EPA  manual  for  wastewater  treatment. 

This  ratio  was  multiplied  by  the  value  for  to  obtain  ~  0.70  d-V  The 

remainder  of  the  coefficients  listed  in  the  table  above  were  (Attained  fiom  a  variety  of 
sources,  including:  (Abbott,  1992),  (Metcalf  and  Etkly,  Inc.,  1991),  and  class  notes 


from  a  course  entitled  "Biological  Treatment  Systems",  H.D.  Stensel,  1993. 
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parameter  definitions,  and  tyfHcal  values  assigned).  For  a  complete  listing  of  mass 
balance  equations,  see  Appendix  C. 


Description  of  Mass  Balance  Equations 


Ammonia  Mass  Balances 

Stage  1  (Ax-1> 


dt 


mw + iQ^s  )iK) + -  (a + Q..S  + 

JTenB  (a)] 


( u  ^ 

f  ) 

y.  Xrs  J 

(K,+N,) 

[T«rai(b)} 


(Tcraitc)] 


(  V,p,„  \ 

yn 

iQoP<,)HQc-^Q^sWo)- 

“J+(a(.)Xi’,) 

[T«ni(4)] 


^K,){X„W.){F,){F^) 


NO, 


NO,-^NO, 

ITenn  (c)] 


Term  (a)  represents  the  flow  of  ammonia  into  and  out  of  the  first  anoxic  zone 
in  the  basin.  Term  (b)  represents  the  depletion  of  ammonia  by  Nitrosomonas 
(nitiifiers).  The  dissolved  oxygen  value  (O])  is  included  to  take  into  account  the 
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inhibiting  effect  of  low  DO  on  nitrification  rates.  In  the  anoxic  stages,  the 
dissolved  oxygen  concentration  (O])  is  assumed  to  equal  zero.  Therefore,  ii.  the 
first  five  stages,  term  (b)  is  equal  to  zero.  Terms  (c)  and  (d)  represent  the 
depletion  of  ammonia  by  the  heterotrophic  uptake  of  soluble  COD  (Sj)  and 
particulate  COD  (P]).  Term  (e)  represents  file  addition  of  ammonia  to  file 
wastev^er  due  to  the  decay  and  deafii  of  heterotrophic  bacteria.  The  value  for 
Fdn  is  included  in  term  (e)  because  only  a  fi’action  of  the  heterotrophic  organisms 
are  able  to  use  nitrate  as  an  electron  acceptor. 


The  ammonia  mass  balance  for  the  first  aerobic  zone  (Ox-2)  is  the  same  as  the 
equation  for  Stage  1  above,  except  that  term  (e)  is  replaced  wifii  the  following: 


HKX.x,W.KfJ 


— 1 
k,+o) 


This  ammonia  production  term  (for  endogenous  decay)  takes  into  consideration  the 
limiting  effects  of  low  dissolved  oxygen  concentration.  The  entire  ammonia  mass 
balance  for  stage  6  (Ox-2)  is  written  as  follows: 


=  0  =  iQo+QKMs 


+  S1kI)'’'Sk3))(^5)”(Co+2#U5  +QkI)'''21(J)K^6) 


[Xen*  (•)! 


rTcni(k)) 

~^h[(Qo  QlUS  ■*"Qk3)K‘^5)~(Co  Qr(3)  X'^is  )](fv  ) 

~fff[(fio  ■^flr<3)K^)”(Co  ■*" i3lr<3) K ■^ )]( i'V  ) 

(Tcm(d)] 
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rr««»(«)l 


=  0  =  (ax^o,)  +  (0^,  )(0.  IXNOJ  +  (Q.r,XNO,)  -  (Q,  +  +  Q.nXNO,) 


term  (a) 


term  (b) 


1.42 

2.86 


*0*  Y  NO, 

k^+0,XNO,-^NO, 


term(c) 


r, 


term(d) 


Term  (a)  represents  the  flow  of  nitrate  into  and  out  of  the  first  anoxic  zone  in 
the  basin.  Term  (b)  represents  the  depletion  of  nitrate  due  to  heterotrophic  uptake 
of  soluble  substrate.  Term  (c)  represents  the  use  of  nitrate  during  endogenous 
decay,  and  term  (d)  represents  nitrate  used  due  to  particulate  COD  degradation. 
Term  (d)  has  been  multiplied  by  fi-actions  that  account  for  flie  inhibiting  effect  of 
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high  DO  or  low  nitrate  concentrations  on  die  rate  of  nitrate  depletion  1^ 
heterotrophic  organisms  using  NO3*  as  an  electron  acceptor. 


Stage  6  (Ox-21 

The  mass  balance  for  nitrate  in  the  first  aeration  zone  (Ox-2)  is  as  follows; 


<mo. 


= 0= (a +0.., +a<„ 


term  (a) 


V  ) 


((X„)(N02.)Y  o,  1. 
K,+N02,  A*-  +  0j‘ 


term(b) 


ri-uiy.Y;/..,,,'!..  (  s,x,  Y  *, 


\F, 

«erm(c) 


2.86  Y„  )^[K.+S.Kk^  +  0. 


Y  NO, 


tenii(d) 


_r V  {Kf)xS  — Y — ^ — ) 
l  2.86  ;  W-^0,)^NO,-^NOj 


tenn(e) 


Term  (b)  is  added  to  this  mass  balance  to  represent  the  production  of  nitrate  from 
nitrite  by  Nitrobacter. 
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Nitrite  Mass  Balances 

Stage  1  (Ax-ll 

= 0 =(ax  w2.)+(a^xA^02,.)+(a,„xAf02,)-(a +a« +a«)(Af02,) 

This  equation  describes  the  flow  of  nitrite  into  and  out  of  stage  1 .  This  mass 
balance  assumes  that  there  is  no  production  or  depletion  of  nitrite  in  the  anoxic 
zones,  and  that  nitrite  is  not  used  for  denitrification  in  the  oxic  zones. 

Stage  6  (Ox-2) 


ydN02, 

dt 


= 0 = (fl, + a., + + a<„)(  W2,)  -  (a + q.^  *  a„> + a<,,K  wo2.) 

term  (a) 


tcnn(b) 


r  X„m2,  ) 

1  y^B  J 

temi(c) 


Term  (b)  represents  the  production  of  nitrite  by  Nitrosomonas.  Term  (c) 
represents  the  depletion  of  nitrite  Nitrobacter. 
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SCOP  Mass  Balances 


Stage  1  (Ax-1) 


=0= as. +e»,A  +a,.,s,  -(a +a„,)(s.)-(fi.)(K) 


>^diere: 


's,x,Y  Y  m 


U.+s, 


A 


k^^O,lNO,^NO, 


Rgu  represents  the  rate  of  depletion  of  the  soluble  COD  (SCOD). 


Stage  6  (Ox-21 


dS. 


Term  (a) 


-^V 


Tenn(b) 


V 


'"ah 


NO. 


k,  +  Sj[k^-^Oj{NO,+NOj 

Term  (c) 


KiF^) 


Term  (a)  represents  the  flow  of  SCOD  into  and  out  of  the  first  oxic  zone  (stage  6). 
Term  (b)  describes  the  depletion  of  SCOD  by  heterotrophic  organisms  using 
oxygen  as  an  electron  acceptor.  A  low  DO  concentration  inhibits  the  SCOD 
utilization  rate.  Term  (c)  represents  the  depletion  rate  of  SCOD  by  heterotrophic 
organisms  using  nitrate  as  an  electron  acceptor.  High  DO  or  low  nitrate 
concontrations  inhibit  the  rate  of  SCOD  depletion. 
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PCOD  Mass  Balances 


Stage  1  (Ax-11 


^^=o=(a)(/>,)+(a+a„.x/i.)-^^i|^+<a„,x^,)-(a+e^+a<„K/i)-/-,(>'.) 


where:  r=F^{KRX^ 


{  Y  NO, 
*-+0,AArO,+AfO,, 


The  term  labeled  describes  the  rate  of  depletion  of  PCOD  by  heterotrophic 
organisms  under  anoxic  conditions. 


Stage  6  (Ox-2t 


o=(a +e,„,  +a,.,  +a«,)(^,)-(a +0u.  +a,., +a,„)(^.) 

Term  (a) 


-{KRX„V,) 


-{KRX„vA 


NO.  '1 

Term(b) 

[no,+noJ 

r  0,  ] 

iF^) 


Term  (c) 


Term  (b)  describes  the  rate  of  depletion  of  particulate  COD  at  low  or  zero  DO 
concentrations  using  nitrate.  Term  (c)  represents  PCOD  degradation  at  higher  DO 
concentiations,  \^ere  oxygen  is  used  as  the  electron  acceptor. 
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Oxygen  Consmn^tion 

Equations  for  the  dissolved  oxygen  consumption  rate  have  been  included  in  the 
five  aerobic  stages  of  this  model  (stages  6  through  10)  in  units  of  kg  O2  /  hr.  A 
present  limitation  of  this  model  is  that  it  does  not  account  for  the  oxygen  used  in 
the  anoxic  stages  due  to  the  use  of  course  and  fine  bubble  aeration  for  mixing. 

The  equation  for  stage  6  (Ox-2)  is  shown  below  as  an  example: 


Stage  6  (Ox-2i 


Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 


(^)(r„)  +  1.42(J^,)| 
term  (a) 

^^ilere: 


O, 


|(F,)+(3.22)(r^„,<„,)(K,)+(UlXr;,^,^„,)(i;) 
term  (b)  term  (c) 


= 


\  j 


+K,(/>)(2r„) 


a 


A  =  l-142(Y„) 


\  ^NB  ) 


xK,+Ng 


O, 


'NOt(ax) 


0,  ] 

<  ^NB  J 

{  K,+N02g  ) 
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Term  (a)  rqjresents  the  oxygen  consumed  for  heterotrophic  degradation  of  soluble 
and  particulate  substrate,  and  for  endogenous  decay.  This  term  also  describes  the 
oxygen  consumed  during  the  endogenous  decay  of  the  heterotrophs.  Tem  (b) 
represents  the  oxygen  consumed  as  a  result  of  depletion  of  ammonia  by 
Nitrosomonas.  Term  (c)  represents  the  oxygen  consumed  due  to  depletion  of 
nitrite  by  nitrobacter  under  aerobic  conditions. 

Approximation  of  Biomass.  Inert  Solids  and  PCOD(]0)  Concentrations 

An  important  and  reasonable  simplifying  assumption  used  to  solve  the  model 
equations  was  that  the  biomass  concentrations  were  constant  from  stage  to  stage. 
This  required  performing  a  steady  state  mass  balance  for  solids  for  the  overall 
system.  A  similar  approach  was  used  for  the  particulate  COD  to  aid  in  the  model 
solution.  Mass  balance  equations  were  done  for  the  entire  ten-stage  system  (as  a 
whole)  to  approximate  the  following  components: 

Xh  (Heterotrophic  biomass  concentration) 

Xns  (Nitrosomonas  Nitrifier  concentration) 

Xnb  (Nitrobacter  Nitrifier  concentration) 

Xj  (Inert  solids  concentration) 

P 1 0  (effluent  PCOD  concentration) 

These  mass  balance  equations  have  been  reazranged,  and  are  listed  below: 

(  rAASXSRT)  ' 
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where: 


AiS  =  5o  ~  ‘S'lo  +  ^ 


[iSRT)a) 


^  _Y,,(N,-NJ(SRT) 


^  _Y,,{N,-NJiSRT) 
""  il  +  k,,JSRT))0 

^  a(X,,,^)(SRT) 


Xt  -  Xh  +  Xus  ^nb  Xp 

Pl0=  (1.3)Xp 

[This  approximation  assumes  the  relationship  between 
PlO  and  Xp  can  be  defined  with  the  ratio: 

CODA/^SS  =  1.3.  This  ratio  is  used  to  convert  Pjo  to  Xp] 

Assume 


P.= 


an 


-^  +  K(XJ(Vr) 

(SRT)  ”  "  ^ 


P(10)  (solve  for  Kp) 


The  value  for  X^  (total  MLSS)  is  obtained  fi’om  Phoenix  Plant  data  for  the  testing 
date  that  is  used  in  the  model.  The  last  equation  can  then  be  rearranged  to  solve 
for  Kp,  with  the  assumption  that  Pio  =  Pavg 
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Incorporatioii  of  Model  into  an  Excel  Snreadshect 


Appendix  A  shows  the  layout  of  the  model  on  the  Excel  spreadsheet. 

Column  (B)  lists  the  calculated  concentrations  for  each  of  the  components  being 
modeled  These  components  are  calculated  using  the  Excel  "solver”  function,  and 
are  computed  for  each  of  the  ten  stages  of  the  anoxic-aerobic  basin.  The  six 
components  included  on  the  spreadsheet  are  as  follows; 


SCOD  — 
NO 

PCOD  — 
N 

N02 

Ox  consu  - 


Soluble  COD  concentration  (mg/1) 
Nitrate-nitrogen  concentration  (mg/l) 
Particulate  COD  concentration  (mg/l) 
Ammonia-nitrogen  concentration  (mg/l) 
Nitrite-nitrogen  concentration  (mg/l) 
Oxygen  consumed  (kg02/hr) 


Columns  (D)  through  (H)  of  the  spreadsheet  contain  the  output  values  produced  by 
mass  balance  equations  for  each  component  in  each  of  the  stages.  These  mass 
balances  have  been  written  in  terms  of  ceil  references.  ^  kinetic  constants  and 
initial  input  values  for  this  model  are  located  in  columns  (L)  through  (X)  of  the 
spreadsheet.  The  lower  portion  of  die  spreadsheet  contains  output  graphs  and 
tables,  describing  the  changing  concentrations  through  the  ten  stage  system  for 
each  component  being  modeled.  Typical  concentration  curves  based  on  plant  data 
are  also  provided  for  comparison  to  the  curves  computed  by  mass  balances. 
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Model  Solution 


Any  of  the  input  values  in  Table  A.  1  of  Appendix  A  [columns  (N)  dirough 
(X)  on  the  Excel  spreadsheet]  may  be  altered  for  tiie  purpose  of  analyzing  the 
effect  on  the  anoxic-aerobic  treatment  system.  When  any  of  these  input  values  are 
changed,  the  "solver  function"  in  Excel  must  be  used  first  to  solve  the  two 
simultaneous  equations  for  and  P(]0)-  These  equations  are  located  in  ro^^-s 
41-42  and  columns  (N)  through  (O)  on  die  spreadsheet.  This  solver  command  can 
be  activated  by  clicking  the  Macro  button  in  cell  0-37. 

The  "solver  command"  must  then  be  used  ten  more  times  for  each  of  the  ten 
stages  in  the  treatment  system.  Using  the  first  anoxic  stage  (Ax-1)  as  an  example, 
the  "solver  function"  would  be  selected  from  the  Format  Menu  of  Excel.  Cell  D6 
would  be  used  as  the  target  cell  in  stage  1.  Cells  D6  through  DIO  would  be  set 
equal  to  zero  as  constraints  in  the  "solver"  dialog  box.  Cells  B6  through  BIO 
would  be  designated  as  the  "changing  cells".  The  table  below  lists  die  input 
requirements  for  the  dialog  box  of  the  "solver  function"  in  each  of  the  ten  stages. 

A  macro  button  has  jieen  installed  for  each  stage  to  speed  up  the  solving  process. 
Clicking  this  button  with  the  mouse  will  automatically  activate  the  "solver 
function"  for  the  corresponding  anoxic  or  aerobic  stage.  The  button  location  for 
each  stage  is  also  listed  on  the  following  table: 
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Table  2;  Input  Requirements  for  Microsoft  Excel 


"Solver"  Function 


Stage 

Target 

Cell 

Changing  Cells 

Constraints 

Macro 

Button 

Location 

1 

D6  =  0 

B6:B10 

D6;D10  =  0 

A5 

2 

D14=0 

B14;B17 

D14;D17=0 

A13 

3 

D22=0 

B22:B26 

D22:D26=0 

A21 

4 

D30=0 

B30:B33 

D30:D33=0 

A29 

5 

D38=0 

B38:B41 

D38:D41=0 

A37 

6 

D46=0 

B46:B50 

D46:D50=0 

A45 

7 

D54=0 

B54:B58 

D54;D58=0 

A53 

8 

D62=0 

B62:B66 

D62:D66=0 

A61 

9 

D70=0 

B70:B74 

D70:D74=0 

A69 

10 

D78=0 

B78:B82 

D78:D82=0 

A77 

There  will  normally  be  at  least  one  stage  within  each  run  in  ^^iiich  the  "solver 
function"  cannot  find  a  feasible  solution  for  the  set  of  simultaneous  equations. 
When  this  occurs,  the  component  that  did  not  achieve  a  successful  solution  must 
be  identified,  and  "re-solved"  individually  using  the  "solver  function”.  This 
component  can  be  identified  by  observing  the  values  listed  in  column  (D).  If  the 
"solver”  has  found  a  successful  solution,  then  all  of  the  values  in  column  (D) 
should  be  very  close  to  zero.  If  any  of  the  values  in  colunm  (D)  are  not  close  to 
zero,  then  the  corresponding  component  must  re-solved.  This  provides  a  seed 
value  that  can  be  used  to  re-solve  all  of  the  mass  balances  for  a  particular  stage 
simultaneously.  The  new  seed  value  [re-entered  into  the  appropriate  row  of 
column  (B)]  should  allow  for  a  successful  solution  using  the  "solver  function"  .  If 
not,  then  repeat  the  procedure  for  unsolved  mass  balances  until  appropriate  seed 
values  are  obtained  (see  Microsoft  Excel  "User's  Guide"  for  description  of 
"Solver"  command). 
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Chapter  3:  Substrate  Storage  Model 

As  previously  stated,  significant  substrate  uptake  and  storage  of  SCOD  occurs 
in  the  first  anoxic  zone  in  the  aeration  basin.  An  additional  model  has  been 
developed  that  accounts  for  substrate  storage  on  an  Excel  SpreadshcK^,  and  is 
referred  to  in  this  paper  as  the  "storage  model".  This  model  accounts  for  the 
initial  accumulation  and  subsequent  utilization  of  stored  substrate  within  the  cell. 
Steady-state  conditions  have  again  been  assumed  for  diis  model.  The 
concentrations  of  the  following  components  are  computed  for  each  of  the  ten 
anoxic-aerobic  stages: 

ratio  of  soluble  stored  substrate  to  heterotrophic  biomass  (S^/Xh) 
anunonia-nitrogen  (N) 
nitrate-nitrogen  (NO) 
nitrite-nitrogen  (N02) 

soluble  substrate  as  chemical  oxygen  demand  (S) 
particulate  substrate  as  chemical  oxygen  demand  (?) 

In  order  to  accoimt  for  cell  storage,  the  value  for  Sg  is  divided  by  Xh  throughout 
this  model  to  create  a  ratio  (Ss/X{{).  This  normalizes  the  stored  substrate 
concentration  by  relating  it  to  the  heterotrophic  biomass  concentration.  The 
following  mass  balances  are  based  on  the  assumption  that  as  the  amount  of  stored 
substrate  per  unit  biomass  increases,  there  will  be  a  higher  degradation  rate. 


25 


Description  of  Mass  Balance  Equations 

The  following  is  a  listing  of  typical  mass  balance  equations  for  the  "substrate 
storage  model".  See  Table  1  for  coefGcient/parameter  definitions  and  typical 
values.  For  a  complete  listing  of  mass  balance  equations,  see  Appendix  C. 


Sc/Xh  Mass  Balances 
Stage  1  (Ax»n 


/ 

+ 

\ 


Term  (a) 


^hJ 


(aS.-fiS,)(0.7)-k,„)(K,) 

Term  (b)  Term  (c) 


where; 


Y  Y  \  JT  \ 


Term  (a)  represents  the  flow  of  stored  substrate  (Sg)  into  and  out  of  the  first 
anoxic  zone.  Term  (b)  describes  the  uptake  of  soluble  COD  into  cell  storage. 
Based  on  observations  at  the  Phoenix  Plant,  this  model  assumes  that  70  %  of  all 
influent  SCOD  is  taken  up  into  storage  in  the  first  anoxic  stage  (Ax-1).  In  the 
mass  balance  equations  for  stages  2  and  3,  this  uptake  percentage  is  assumed  to  be 


10  %  and  100%  respectively.  Residual  non-biodegradable  SCOD  is  subtracted 
from  the  influent  total  SCOD  concentration.  Term  (c)  represents  the  rate  of 
depletion  of  stored  SCOD. 

The  stored  substrate  mass  balance  equations  for  the  aerobic  stages  (zones  6 
through  10)  are  similar  to  the  above  equation,  with  the  exception  that  the  RssU 
equation  is  substituted  with  the  following: 


One  potential  weakness  in  this  model  is  die  lack  of  information  available  for 
estimating  the  following  coefficients: 

Ksto  (coefficient  for  stored  substrate  utilization) 

kss  (half-saturation  constant  for  stored  substrate  utilization) 

For  the  run  of  the  spreadsheet  model  shown  in  Appendix  B,  values  of  0.937  1/day 
and  0.002  g/g  were  used  for  K^to  and  1^  re^ctively.  However,  the  accuracy  of 
these  values  is  not  known.  The  value  for  K^to  estimated  using  an  iterative 
process  with  the  Ss/X}|  mass  balance  equations  in  stages  5  through  10.  A  K^to 
value  was  chosen  that  would  satisfy  the  condition  that  in  stage  10  is  less  than  or 
equal  to  one  one-hundredth  (1/100)  of  in  stage  5.  In  equation  form: 

Ss(10)  <or=  (0.01)(Ss(5)) 

This  was  based  on  the  assumption  that  the  rate  of  stored  substrate  utilization 
would  rtqhdly  start  to  decrease  in  the  re^on  of  the  curve  Where  Sg^jo)  is  less  than 
(0.01)(Ss(s)).  Oxygen  consumption  data  is  expected  to  be  available  from  the 
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Phoenix  Plant  in  1994,  and  this  will  hopefully  allow  for  a  better  ^}proximation  of 
"the  Ksto  value. 


SCOP  Mass  Balances 

The  SCOD  mass  balance  for  stage  1  in  the  "Storage  Model"  is  shown  below. 
Stage  1 


=  0  =  )(S^^)  -(Qo  +  Qkd  Qms  )(*5|) 

Tcnii  (a) 


-(a5.-as,)(o.70) 

Tcmi(b) 


Term  (b)  assumes  70  %  uptake  into  cell  storage  of  soluble  degradable  COD. 
For  the  mass  balance  equations  in  stages  2  and  3,  10%  and  100%  uptake  is 
assumed  respectively.  Therefore,  after  the  diird  anoxic  stage,  all  degradable 
SCOD  will  theoretically  have  been  dej^eted.  These  SCOD  mass  balances  assume 
that  none  of  the  soluble  substrate  is  metabolized  by  the  microorganisms.  All 
SCOD  that  is  removed  from  die  wastewater  in  the  first  three  stages  is  taken  up  into 
cell  storage.  After  the  third  anoxic  zone  (Ax-3),  all  that  remains  is  the  residual 
non-biodegradable  SCOD  (S^,  which  remains  at  a  constant  concentration  through 
the  remaining  seven  stages  in  the  basin. 


The  remainder  of  the  storage  model  mass  balances  are  identical  to  the 
equations  listed  for  the  model  in  Chsqiter  2  (for  NH4,  NO3,  N02‘,  PCOD,  and 
oxygen  consumption),  with  the  exception  of  the  following  three  modifications: 


28 


Modification  #1 

For  NH4  mass  bailees,  term  (c)  is  replaced  with  the  following  (uang  stage  1 
and  6  as  mi  example): 

i^iiere: 


^5.V  ~  ^0 


„  Y  m  U  . 


{in  anoxic  zones} 


and 


^,V  ~  ^tto 


\r 


yKHS^Jx„)j 


O. 


K-^0, 


{in  aerobic  zones} 


This  portion  of  the  ammonia  mass  balance  represents  die  depletion  of  NH4-N  used 
for  cell  synthesis  during  heterotrophic  utilization  of  stored  soluble  substrate. 


Modification  #2 

For  NO3  mass  balances,  term  (b)  is  replaced  with  the  following  (using  stage  1 
as  an  example): 


where; 


^,V  ~  ^mo 


SmIX,  1 

f  Y  NO,  ^ 

Kk^-^0,lNO,^NO,} 

Fos) 
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The  basic  forni  of  this  term  can  be  used  in  the  NO3  mass  balance  equations  for 
both  the  anoxic  and  aerobic  zones  (all  ten  stages).  This  portion  of  the  ammonia 
mass  balance  represents  the  depletion  of  NO3-N  by  heterotrophic  utilization  of 
stored  soluble  substrate. 

Modification  #3 

For  dissolved  oxygen  consumption,  the  heterotrophic  substrate  utilization  rate 
(rsu)  is  replaced  with  the  following  equation  for  stored  substrate  utilization  rate: 


V 


a 


+  k,(P,)(X„) 


O, 


k^  +  Oj 


Incornoratilon  of  Model  into  Excel  Spreadsheet 

Appendix  B  shows  die  layout  of  die  storage  model  on  the  Excel  spreadsheet. 
An  additional  component  (Ss/X^)  has  been  added  under  column  (B)  of  each 
stage,  along  with  the  associated  output  value  produced  by  mass  balance  equations 
in  columns  (D)  through  (F).  A  chart  displaying  die  fluctuation  of  the  Sg/X^  value 
in  each  of  the  ten  stages  has  also  been  added  on  the  lower  portion  of  the 
spreadsheet.  The  remainder  of  the  spreadsheet  layout  is  identical  to  the  model 
previously  described  in  Ch^qiter  2. 
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Model  Solution 


The  procedures  for  operation  of  this  storage  model  are  similar  to  those 
described  in  Ch^ter  2  of  this  paper.  The  only  difference  is  that  one  additional 
mass  balance  equation  (for  S^/Xh)  must  be  incorporated  into  the  solution  process 
for  each  stage.  Table  3  lists  the  input  requirements  for  the  dialog  box  of  the 
"solver  function"  in  each  of  the  ten  stages  (for  manual  use  of  the  "solver" 
command).  Macro  buttons  have  again  been  installed  for  each  stage  to  speed  up  the 
solution  process. 


Table  3:  Input  Requirements  for  Microsoft  Excel 


”Solver**  Function  (for  Storage  Model) 


Stage 

Target 

Cell 

Changing  Cells 

Constraints 

Macro  Button 
Location 

1 

D6=0 

B6:B10;B12 

D6:D10;D12=0 

A5 

2 

D14=0 

B14:B17;B20 

D14:D17;D20=0 

A13 

3 

D22=0 

B22:B26;B28 

D22:D26;D28=0 

A21 

4 

D31=0 

B31:B33;B36 

D31:D33;D36=0 

A29 

5 

D39=0 

B39:B41;B44 

D39:D41;D44=0 

A37 

6 

D47=0 

B47;B50;B52 

D47:D50;D52=0 

A45 

7 

D55=0 

B55:B58;B60 

D55:D58;D60=0 

A53 

8 

D63=0 

B63:B66;B68 

D63:D66;D68=0 

A61 

9 

D71=0 

B71:B74;B76 

D71:D74;D76=0 

A69 

10 

D79=0 

B79:B82;B84 

D79:D82;D84=0 

A77 
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Chapter  4:  Results  and  Discussion 


Introductiop 


The  storage  model  was  run  using  the  typical  kinetic  coefficients  shown  in 
Table  1 .  The  spreadsheet  results  are  provided  in  Appendix  B.  Profile  data  fi-om 
the  Phoenix  Treatment  Plant  available  for  eight  different  testing  dates  between 
October  1992  to  April  1993.  A  data  set  fi’om  17  November  1992  v;as  selected  for 
analysis  in  this  paper  because  the  concentrations  of  the  various  components  were 
in  the  average  range  v^en  compared  to  the  other  data  sets. 

Some  coefficient  and  parameter  values  fi-om  the  Phoenix  Plant  profile  data 
were  input  directly  into  the  model.  These  include  values  for: 

Qo  (influent  flow)  =  20  MOD 
Qras  (return  flow)  =  9.34  MOD 
Xj  (total  MLSS)  =  3910  mg/1 
Xjis  (nitrosomonas  concentration)  =  62  mg/1 
Xnb  (nitrobacter  concentration)  =  31  mg/1 

The  values  for  ^ns  and  X||b  were  computed  in  a  summary  report  on  the  Phoenix 
Plant  operations(Stensel,  1 993(c)).  These  calculations  were  based  on  the  average 
influent  TKN,  SRT,  and  COD  removed  over  a  two  week  period.  The  original 
intention  was  to  use  the  following  equations  for  calculation  of  Xq^  and  Xq^: 

^  Y,s(N,-N,,-N^}(SRT) 

u+k,,„,(SRT)]e 
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y.BiK-N^o-N^yiSRT) 

[l+k,^„,(SRT)]$ 


where: 


^J[rEr)_^0  09) 
«"  (SRTX&y  ’ 


*** 


Nsyn  =  amount  of  influent  nitrogen 
used  for  cell  growth 


However,  use  of  these  equations  yielded  values  that  were  15  percent  higher  than 
the  average  values  obtained  directly  from  available  plant  data.  Therefore, 
concentrations  of  nitrosomonas  and  nitrobacter  (X^s  and  were  fixed  for  the 
purpose  of  this  analysis  (for  the  SRT  analysis,  the  Xn^  and  Xni)  values  were 
prorated  to  account  for  changing  SRT). 

Table  1  shows  that  there  are  many  different  coefficients  that  could 
conceivably  be  adjusted  to  achieve  a  "better  fitting  curve".  However,  for  the 
purposes  of  this  analysis,  the  only  coefficients  that  were  considered  for  adjustment 
in  curve  fitting  are  as  follows: 

(nitrogen  half  saturation  constant  for  nitrobacter)  =0.6  mg/1 
(maximum  specific  growth  rate  for  nitrobacter  nitrifiers)  =  0.7  d"^ 

Xjo  (influent  inert  solids  concentration)  =  20  mg/1 

The  values  for  KN(ns)  were  fixed  at  0.5  mgA  and  0.5  d"l  based  on 

testing  data  obtained  from  Phoenix  Plant  (Stensel(b),  1993). 
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In  order  to  run  the  spreadsheet  model,  several  initial  assumptions  were  made 
for  influent  concentrations,  and  concentrations  in  stages  9  and  10.  The  internal 
recycle  withdrawal  is  located  in  stage  9,  and  therefore  required  an  assumed 
concentration  for  each  component.  Assumed  concentrations  for  stage  10  were 
needed  due  to  the  return  flow  that  is  cycled  through  a  secondary  clarifier. 

Assumed  concentrations  for  the  following  components  are  taken  directly  fi'om  the 
Phoenix  Plant  Data  (for  the  particular  testing  date  being  analyzed): 

SCOD  (influent) 

SCOD  (stage  9) 

SCOD  (stage  10) 

Nitrate  (stage  10) 

Nitrate  (stage  9) 

Ammonia  (influent) 

Ammonia  (stage  10) 

Ammonia  (stage  9) 

Nitrite  (influent) 

Nitrite  (stage  10) 

Nitrite  (stage  9) 

These  values  were  used  to  initiate  the  computations  in  the  spreadsheet.  The 
calculated  values  (computed  in  the  model)  were  normally  very  close  to  the 
concentrations  that  were  initially  assumed.  Influent  nitrate  was  always  assumed  to 
equal  zero  in  this  model,  and  the  influent  PCOD  concentration  was  assumed  to 
equal  the  influent  SCOD  concentration  (as  a  rough  approximation). 

Figui  es  2  through  8  show  the  results  for  the  storage  model  using  coefficient 
and  parameter  values  as  listed  in  Appendix  B.  The  calculated  NH4,  NO3,  and 
N02“  concentrations  (produced  by  the  model)  correspond  fairly  well  to  the 
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Phoenix  Plant  ProjQle  Data  in  each  stage  of  the  ten  stage  system.  The  calculated 
concentration  for  NH4  decreases  steadily  in  the  aerobic  zones  (stages  6  through 
10)  as  NH4  is  transformed  into  N02"  and  NO  3  through  the  nitrification  process. 
As  shown  in  Figure  2  for  ammonia,  the  calculated  influent  concentration  is  25  % 
greater  than  the  influent  Plant  Data  concentration.  This  difference  is  due  to  the 
assumption  that  the  influent  flow  for  the  Storage  Model  profile  for  ammonia 
includes  TKN  (Total  Kjeldahl  Nitrogen).  However,  the  Phoenix  Plant  Data  only 
includes  the  inorganic  NH4-N,  as  evidenced  by  historical  Plant  Data.  The  influent 
TKN  value  was  assumed  to  be  approximately  25%  greater  than  the  inorganic 
NH4-N  influent  value.  M  of  the  organic  nitrogen  is  assumed  to  be  hydrolyzed  in 
the  first  anoxic  zone  (stage  1).  Consequently,  the  calculated  (model)  curve  and  the 
Plant  Data  curve  are  directly  comparable  in  later  stages  (Figure  2),  due  to  the  fact 
that  most  of  the  inorganic  NH4-N  should  be  released  as  the  number  of  stages 
increases. 

The  calculated  curve  for  Nitrite  (Figure  4)  does  not  correspond  to  the  Plant 
Data  curve  as  well  as  for  NH4  and  NO3.  The  Plant  Data  profile  shows  nitrite 
depletion  in  stages  1  and  2.  The  mass  balance  equations  for  the  Storage  Model  do 
not  include  depletion  terms  for  N02‘  in  the  anoxic  stages.  Therefore,  the  model 
curve  rises  well  above  the  Plant  Data  curve  in  stages  1  and  2.  The  remainder  of 
the  stages  correspond  fairly  well. 

The  model  and  plant  data  curves  for  SCOD  (Figure  5)  fit  very  well,  with  the 
assumption  that  70%  of  the  degradable  SCOD  is  taken  up  into  cell  storage  in  the 
first  anoxic  zone.  The  figure  shows  that  most  of  the  biodegradable  SCOD  is 
depleted  after  reaching  the  third  anoxic  zone,  leaving  only  the  non-biodegradable 
(residual  SCOD)  in  the  remaining  seven  stages. 
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There  is  no  plant  data  available  for  comparison  wiien  evaluating  the  profile 
curves  for  PCOD,  oxygen  consumption  rate,  and  Sg/Xj!  (ratio  of  stored  SCOD  to 
heterotrohic  biomass  concentration).  The  PCOD  curve  is  shown  in  Figure  6.  The 
PCOD  concentration  accumulates  in  stage  1  (due  mostly  to  the  internal  recycle 
flow  from  stage  9),  and  is  then  gradually  degraded  in  stages  2  through  10. 

Figure  7  is  a  bar  chart  showing  the  oxygen  consumption  rates  computed  by  the 
model  in  stages  6  throu^  10.  The  rates  are  computed  in  units  of  kg/hr,  and 
converted  to  mg/l-hr  for  plotting  on  the  bar  chart.  Figure  7  shows  that  the  O2 
consumption  rate  is  fairly  constant  in  stages  6  through  9,  and  drops  slighfly  in 
stage  10  due  to  declining  NH4  and  N02"  concentrations.  No  oxygen 
consumption  is  assumed  to  occur  in  the  anoxic  stages  (zones  1  dirough  5). 

The  stored  substrate  curve  (Sg/Xn)  is  shown  in  Figure  8.  Alfliough  there  is  no 
plant  data  available  for  comparison,  the  curve  behaves  as  would  be  expected.  The 
peak  amount  of  SCOD  in  cell  storage  occurs  in  the  first  anoxic  zone,  with  gradual 
depletion  to  near  zero  in  st^e  ten. 


Analysis  of  Parameter  Adjustments  (Using  Storage  Model) 

The  Storage  Model  as  shown  in  Appendix  B  was  analyzed  to  determine  the 
effect  of  adjusting  various  operating  parameters.  The  following  adjustments  were 
made  (individually)  using  the  coefficients  and  parameters  listed  in  Appendix  B  as 
the  default  values: 

1)  Vary  SRT  (Solids  Retention  Time) 

2)  Vary  DO  (dissolved  oxygen  concentration)  in  aerobic  stages 
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(6  through  10) 

3)  Change  the  internal  recycle  withdrawal  location  from  stage  9  to  stage  10 

4)  Vary  the  percent  distribution  of  internal  recycle  flow  from  stage  9  to 
stages  lands. 

5)  Increase  the  internal  recycle  flow  from  stage  9 

Effect  of  SRT 


The  effect  of  SRT  was  varied  from  3,  5,  7  and  9  days.  The  effect  of  raising 
the  SRT  on  NH4,  NO3  and  N02‘  concentration  can  be  seen  in  Figures  9, 10,  and 
11.  As  SRT  increases,  the  effluent  NH4  concentration  decreases  and  die  effluent 
NO3  concentration  increases  (see  Table  4).  This  can  be  attributed  to  die  fact  that 
the  nitrifier  biomass  increases  as  SRT  increases,  thaeby  resulting  in  a  greater 
nitrification  rate,  and  more  rqiid  depletion  of  NH4.  Figure  1 1  shows  that  the 
effluent  NO2"  concentration  drops  sharply  as  SRT  increases,  due  to  decreased 
amounts  of  NH4  available  for  transformation  to  N02"  by  nitrosomonas  in  stage 
10..  The  SCOD  concentration  remains  fairiy  constant  as  SRT  is  increased,  while 
peak  stored  substrate  (Sj/Xh)  ratios  decline  (see  Figure  12).  The  latter  effect  is 
mainly  due  to  the  increases  in  heterotrophic  biomass  concentration  that  occurs  as 
SRT  increases. 

Effect  of  Dissolved  Oxygen  Concentnition 


The  dissolved  oxygen  concentration  was  varied  in  the  aerobic  zones  (stages  6 
through  10).  DO  concentrations  of  0.5, 1.0, 2.0,  3.0  mgd  were  used.  The  effect 
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of  raising  DO  levels  on  NH4  and  NO3  concentration  profiles  can  be  seen  in 
Figures  13  and  14.  As  DO  is  decreased,  effluent  NH4  concentrations  increase  and 
effluent  NO3  concentrations  decrease  (see  Table  4).  This  change  is  to  be  expected 
due  to  the  inhibitory  effects  of  low  dissolved  oxygen  concentrations  on  the 
nitrification  rate.  As  the  DO  level  falls  below  2  mg/1,  the  rate  in  v^diich  NH4  is 
converted  into  N02"  and  NO3  is  reduced,  resulting  in  higher  effluent  NH4 
concentrations.  There  is  little  change  in  the  profile  for  nitrite,  with  only  a  small 
decrease  in  effluent  NO2'  concentration  in  stage  10  as  DO  increases. 


Effect  of  Internal  Recycle  \^thdrawal  Location 

With  the  current  aeration  basin  design,  the  internal  recycle  flow  is  withdrawn 
from  stage  9.  One-third  of  this  flow  is  directed  to  stage  1  (20  MOD),  and  the 
remainder  (40  MOD)  feeds  into  stage  3.  As  part  of  this  analysis,  the  Storage 
Model  was  modified  so  that  the  internal  recycle  flow  was  withdrawn  fi’om  stage  10 
instead  of  stage  9.  The  effect  on  the  NH4,  NO3  and  N02"  concentration  profiles 
can  be  seen  in  Figures  IS  through  17.  Changing  to  stage  10  withdrawal  lowers  the 
NH4  and  N02"  curves  slightly.  This  is  due  to  the  lower  concentrations  of 
ammonia  and  nitrite  that  are  being  recycled  fium  stage  10  (in  comparison  to  stage 
9).  The  change  in  withdrawal  location  raises  the  NO3  profile  (Figure  16)  due  to 
the  higher  nitrate  concentration  that  exists  in  stage  1 0.  However,  the  effluent 
concentrations  of  each  component  remain  approximately  the  same  (see  Table  4). 
This  would  suggest  that  little  benefit  would  be  realized  as  a  result  of  relocating  the 
withdrawal  location.  There  is  little  change  observed  in  the  SCOD,  PCOD  and 
stored  substrate  curves  as  a  result  of  changing  the  withdrawal  location. 
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Effect  of  Internal  Recycle  Flow  Distribution 

The  internal  recycle  flow  distribution  to  stages  1  and  3  (from  stage  9)  was 
varied  with  five  different  scenarios  as  follows: 


Qir(l) 

Qir(3) 

60  MGD  =  20  MGD 

40  MGD 

30  MGD 

30  MGD 

40  MGD 

20  MGD 

50  MGD 

10  MGD 

60  MGD 

0  MGD 

As  can  be  seen  in  Figures  18  and  19,  the  re-distribution  of  internal  recycle  flow 
does  not  have  a  significant  effect  on  the  NII4  and  NO3  profiles,  except  for  in 
stages  1  and  2.  As  the  percentage  of  internal  recycle  flow  into  stage  1  is  increased, 
the  concentration  of  NH4  in  stages  1  and  2  is  decreased,  and  the  concentrations  of 
NO3  and  N02~  in  stages  1  and  2  are  increased.  The  efiQuent  concentrations  for  all 
five  distribution  scenarios  remain  fairly  constant  (see  Table  4).  SCOD 
concentrations  drop  slightly  in  stages  1  and  2  as  increased  flow  is  fed  into  stage  1 
(SCOD  declines  by  approximately  5%  as  flow  into  stage  1  is  increased  from  20 
MGD  to  60  MGD).  Similarly,  PCOD  concentration  also  fell  slightly  in  the  first 
two  anoxic  stages,  with  the  effluent  levels  remaining  essentially  die  same. 

Figure  20  shows  the  effect  of  varying  internal  recycle  distribution  on  stored  SCOD 
concentration,  with  declining  S^/Xh  values  as  flow  to  stage  1  is  increased. 
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Effect  of  Increased  Internal  Recycle  Flow 

The  existing  internal  recycle  flow  from  stage  9  is  ^jproximately  60  MGD. 

The  Storage  Model  was  used  to  analyze  the  effects  of  varying  this  withdrawal  flow 
as  follows: 


Qir(l) 

Qir(3) 

40  MGD  = 

13.33  MGD 

26.67  MGD 

60  MGD  = 

20  MGD 

40  MGD 

75  MGD  = 

25  MGD 

50  MGD 

Each  of  these  three  scenarios  assumes  that  one-third  of  the  internal  recycle  flow  is 
directed  into  stage  1 ,  with  the  remainder  being  fed  into  stage  3.  Figures  21  and  22 
show  the  effect  of  increased  flow  on  the  ammonia  and  nitrate  concentrations.  As 
internal  recycle  flow  increases,  the  NH4  concentration  decreases  in  stages  1 
through  7,  while  NO3  concentrations  increase  in  stages  1  through  7,  and  decrease 
through  the  remainder  of  the  aerobic  zones.  Table  4  shows  that  NH4  effluent 
concentrations  change  very  little  with  increased  flow.  Little  change  occurs  in  the 
N02~  profile,  with  only  small  concentration  increases  in  die  anoxic  zones  as 
internal  recycle  flow  increases.  The  stored  SCOD  curve  lowers  only  idightly  as 
flow  increases. 


40 


Analvsis  of  Alternate  Data  Set  (Phoenix  Plant  Data  from  4  November  1^2 


The  spreadsheet  run  shown  in  Appendix  B  was  based  on  Phoenix  plant  data 
dated  1 7  November  1 992.  In  order  to  determine  how  well  the  Storage  Model 
would  "fit”  to  another  set  of  plant  data,  an  additional  nm  was  performed  using  the 
4  November  1992  data  set  The  results  are  shown  in  Figures  23  through  29.  As 
can  be  seen,  the  ammonia  and  nitrate  profiles  do  not  fit  as  well  as  was  observed 
widi  the  17  November  Plant  Data.  However,  there  still  seems  to  be  an  adequate 
coirelation  with  the  Plant  Data  curve.  The  SCOD  model  curve  (Figure  26)  falls 
below  the  Plant  Data  curve,  suggesting  that  die  percentage  of  SCOD  uptake  into 
cell  storage  is  significantly  less  than  the  assumed  70%  (in  stage  1).  Figure  30  is 
{Hovided  to  show  the  wide  variability  in  NH4  Plant  Profile  Data  for  different 


ComDarison  of  Michaelb-Menton  Model  to  Storage  Model 


The  Michaelis-Menton  Model  (described  in  Chq)ter  2)  has  been  run  using  the 
same  coefficients  and  parameters  used  for  die  Storage  Model  (Appendix  B).  The 
results  of  the  Michaelis-Menton  Model  are  provided  in  Appendix  A,  including 
concentration  profile  charts  for  each  component  being  modeled.  The  curves  for 
NH4  and  N02~  are  similar  to  the  Storage  Model  in  terms  of  coire^xmdence  widi 
the  Plant  profile  data.  The  PCOD  concentration  profiles  for  the  two  models  are 
also  very  similar.  The  NO3  profile  in  the  Nfichaclis-Menton  model,  however,  is 
considerably  lower  than  for  the  Storage  Model  (see  Figure  31).  The  calculated 
SCOD  curve  in  the  Michaelis-Menton  Model  (Figure  32)  does  not  take  into 
account  non-biodegradable  SCOD.  This  explains  the  poor  correspondence 
between  the  model  and  Plant  Data  curves  for  SCOD.  This  can  be  compared  with 
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the  SCOD  profile  for  file  Storage  Model  (Figure  5),  which  does  take  the  residual 
(non-biodegradable)  SCOD  into  account,  and  provides  a  much  closer 
approximation  of  the  actual  concentrations  in  each  stage.  The  SDNR  (Specific 
Denitrification  Rate)  curve  shown  for  the  Michaelis-Menton  Model  in  Ai^)endix  B 
fits  much  better  to  the  Plant  Data  than  die  SDNR  curve  for  the  Storage  Mottel 
(shown  in  Appendix  A).  This  suggests  that  the  Storage  Model  is  not  suffidendy 
calibrated  due  to  lack  of  information  for  kinetic  coeffidents. 
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Chapters:  Conclusions 

Two  models  have  been  devdoped  to  emulate  various  components  of  an 
anoxic-aerobic  treatment  system  in  Phoenix,  Arizona.  An  initial  observation 
suggests  that  both  models  sqipear  to  be  equally  suffident  in  thdr  ability  to 
^jproximate  NH4,  N03~,  and  N02‘  concentrations  in  each  of  the  ten  sta^s. 

When  comparing  modeled  concmtrations  to  Phoenix  Plant  Data,  neither  model 
produces  a  predse  estimate  of  actual  concentrations  widiin  the  aeration  basin. 

Only  a  rough  approximation  is  obtained.  The  Storage  Model  should  provide  a 
more  reasonable  correlation,  due  to  the  fact  that  it  takes  into  consideration  the 
uptake  of  soluble  degradable  COD  into  cell  storage.  The  Nhchaelis-Menton 
Model  does  not  take  this  into  account.  However,  in  view  of  the  results,  no 
conclusion  can  be  reached  concerning  die  superiority  of  one  model  over  the  other. 

An  analysis  was  performed  using  the  Storage  Model  to  determine  the  effects 
of  adjusting  various  operating  parameters.  The  model  responded  in  a  logical 
maimer,  predicting  decreased  NH4  concentrations  and  increased  NO3 
concentrations  as  the  SRT  and  DO  levels  were  increased.  The  internal  rer^cle 
withdrawal  location  was  changed  from  stage  9  to  stage  10  as  part  of  the  analysis. 
The  model  showed  that  this  modification  would  have  little  effect  on  die  efQuent 
NH4,  NO3,  and  N02'‘  concentrations.  However,  the  model  could  not  be  calibrated 
due  to  lack  of  certainty  in  the  kinetic  coefiBcients.  Appropriate  values  for  the  K^to 
and  kgs  storage  terms  are  not  known,  and  had  to  be  roughly  approximated  for  the 
purposes  of  this  model.  Therefore,  the  model  caimot  be  used  to  predict  profile 
concentrations  of  the  various  components.  Nevertheless,  die  Storage  Model  could 
be  used  to  predict  trends  in  profile  concentrations  as  a  result  of  changing  various 
oper^ng  parameters. 
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Ammonia  Concentration 


FIGURE  2:  Storage  Model;  Ammonia  Profile 
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Stage 
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FIGURE  3:  Storage  Model ;  Nitrate  Profile 
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FIGURE  4:  Storage  Model;  Nitrite  Profile 
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FIGURE  6:  Storage  Model;  PCOD  Profile 
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Oxygen  Consumption 
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Stage 


FIGURE  7:  Storage  Model;  Oxygen  Consumption  Rate 


Ss/Xh  (mg/L/mg/L) 


Stored  Substrate  Concentration 


FIGURE  8:  Storage  Model;  Stored  SCOD  Profile 
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TABLE  4:  Effluent  Concentrations  of  NH4,  NO3  and  N02‘ 
EFFLUENT  CONCENTRATIONS 
(Storage  Model) 


Ammonia 

Nitrate 

Nitrite 

(mg/1) 

(mg/1) 

(mg/1) 

Internal  Recycle 

Withdrawal  Location 

Stage  9 

2.32638 

10.5886 

0.95286 

Stage  10 

2.7225 

9.78687 

0.98328 

SRT 

3  days 

5.98667 

7.70259 

1.10224 

5  days 

2.32638 

10.5886 

0.95286 

7  days 

0.56086 

12.6163 

0.48776 

9  days 

0.16648 

13.6194 

0.199 

DO 

0.5  mg/1 

6.08655 

5.66895 

1.0997 

1.0  mg/1 

4.21585 

8.12662 

1.07492 

2.0  mg/1 

2.85036 

9.89541 

1.00315 

3.0  mg/1 

2.32638 

10.5886 

0.95286 

Dbtribution  of 

Internal  Recycle  Flow 

Qir(ll .  Oir(3) 

20  MOD  40  MOD 

2.32638 

10.5886 

0.95286 

30  MOD  30  MOD 

2.32554 

10.5291 

0.95274 

40  MOD  20  MOD 

2.32546 

10.507 

0.95273 

50  MOD  10  MOD 

2.32557 

10.4987 

0.95274 

60  MOD  0  MOD 

2.32573 

10.4966 

0.95276 

Increased  Internal 

Recycle  Flow 

40  MOD  =  13.33+26.67 

2.27488 

11.3958 

0.95736 

60  MOD  =  20  +40 

2.32638 

10.5886 

0.95286 

75  MOD  =  25  +50 

2.35185 

10.3736 

0.94915 

51 


TABLE  4  (continued):  Effluent  Concentrations 


Ammonia 

Nitrate 

Nitrite 

Comnare  Storage  Model 
and  Michaelis-Menton 
Model 

Storage  Model 

2.32638 

10.5886 

0.95286 

Michaelis-Menton  Model 

2.39589 

9.23093 

0.95974 

Comnare  4  NOV  and 

17  NOV  Plant  Data  Using 
Storage  Model 

17  NOV 

2.32638 

10.5886 

0.95286 

4  NOV 

2.72203 

11.8047 

0.97504 

52 


cone,  of  NOS  (mg/I) 


53 


Nitrate  Concentration 


- • — 

-  SRT=ad 

- <r— 

-  SRT=6d 

. ♦ 

-  SRT=7d 

- B — 

-  SRT=8d 

Stage 


FIGURE  10:  Effect  of  SRT  on  Nitrate  Concentrations 


(mg/I) 


Nitrite  Concentration 


FIGURE  11:  Effect  of  SRT  on  Nitrite  Concentrations 


Stored  Substrate  Concentration 


FIGURE  12:  Effect  of  SRT  on  Stored  SCOD  /  Xh  Ratio 


56 


FIGURE  13:  Effect  of  DO  on  Ammonia  Concentration 


cone,  of  NOS  (mg/I) 


Nitrate  Concentration 


FIGURE  14:  Effect  of  DO  on  Nitrate  Concentration 


cone,  of  NH4  (mg/I)) 


38 


FIGURE  15:  Effect  of  Internal  Re^rcle  Withdrawal  Location 
on  Ammonia  Concentration 


FIGURE  16:  Effect  of  Internal  Recycle  Withdrawal  Location 
on  Nitrate  Concentration 


60 


Nitrite  Concentration 


« —  Stage  9Witidrav«l  — a —  Sttge  10Wiihdrav«l 


FIGURE  17:  Effect  of  Internal  Recycle  Withdrawal  Location 
on  Nitrite  Concentration 


cone,  of  NH4  (mg/I) 
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FIGUIIE  18:  Effect  of  Internal  Recycle  Flow  Distribution 


on  Ammonia  Concentration  [  Qir(l)  /  Qir(3)  ] 


cone,  of  NOS  (mg/I) 


62 


FIGURE  19:  Effect  of  Interaal  Recycle  Flow  Distribution 


Nitrate  Concentration  [  Qir(l)  /  Qir(3)  ] 


FIGURE  20:  Effect  of  Internal  Recycle  Flow  Distribution 

on  Stored  SCOD  Profile  [  Qir(l)  ^  Qir(3)  ] 


cone,  of  NH4  (mg/I) 
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FIGURE  21:  Effect  of  Internal  Recycle  Flow  Amount 


on  Ammonia  Concentration 
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cone,  (mg/l) 


Ammonia  Concentration 


FIGURE  23:  4  November  1992  Data  Set;  Ammonia  Concentration; 

Storage  Model 
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Storage  Model 


FIGURE  25:  4  November  15^2  Data  Set;  Nitrite  Concentration; 

Storage  Model 
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SCOO  Concentration 


Stage 


« —  CalculaiBd  — * —  PlantData 


FIGURE  26:  4  November  1992  Data  Set;  SCOD  Concentration; 

Storage  Model 


FIGURE  27:  Compare  4  November  and  17  November  1992  Data  Sets; 
PCOD  Concentration;  Storage  Model 
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Oxygen  Consumption 


60^ 


1  23456789  10 

Stage 


FIGURE  28:  4  November  1992  Data  Set;  Oxygen  Consumption  Rate; 

Storage  Model 


FIGURE  29:  4  November  1992  Data  Set;  Stored  SCOD  ; 


Storage  Model 


FIGURE  30:  Variation  of  Anunonia  Concentration  Profiles 

from  Different  Phoenix  Plant  Data  Sets 


FIGURE  31:  Michaelis-Menton  Model;  Nitrate  Concentratioii 


conc.(mg/l) 


SCOD  Concentration 


FIGURE  32:  Michaelis-Menton  Model;  SCOD  Concentration 
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Appendix  A:  Model  Based  on  Michaelis-Menton  Kinetics 


78 


79 


80 


E  F  G  H 


(MICHAELIS  -  MENTON  KINETICS  MOD 


SCOD  = 


\\M 


8  PCOD  = 


9  N 


10  N02  = 


11  Ox  consu 


12 


13 


14  SCOD  = 


15  NO 


16  PCOD  = 


17  IN 


19  Ox  consu 


AX-1 


68.06068 


480.06 


20.89265 


0.537893 


AX-2 


65.65949 


0.073741 


479.8755 


20.8045 


-1.3E-07  3853.882 


-2.4E-08  105.7365  2.896626|  8.060021  I 


-5.3E-09  293.7714 


-1E-06i  01  21.0422]  1.603992|  0.583405| 


3.55E-14 


2.13E-07 


3.79E-08  20.84084 


2.33E-09  62.78873 


4.147451 

0.318653 

0.1 159451 

21  1  \AX- 


SCOD  =  23.21862 


23  NO  =  ;  0.423677 


24 


25  IN  =  i  13.82667 


27  Ox  Consu 


-2.  IE-08 


-IE-06 


-IE-06 


205.0185  8.079959  20.68507 


0  40.79991  4.116452  1.627372 


SCOD  = 


31  |NO  = 


33 


34  N02  = 


35  I  Ox  Consu 


SCOD  = 


0 

PCOD  = 


imu. 


AX-4 


I  21.82914 


0.144015 


441.5274 


□I 


AX-S 


i  21.3406 


441.4762 


02  = 


Ox  Consu 


-7.1E-081  1001.212 


21.834821  0.8858771  2.2671461 


-2.9E-07I  103.9575 


-IE-06 


4.3452591 

0.451175 

0.178423! 

-1.7E-08  352.0229  I 


7.6770531  0.314881  0.805753 


-1.9E-07  36.94691 


-IE-06  0  1.527779  0.16035  0.06342! 
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A 

B 

C 

0 

E 

F 

G 

H 

1 

o 

m 

EH 

SCOD  = 

4.657131 

-7.9E-11 

1477.689 

12.97904 

NO  = 

0.835624 

-3.4E-09 

73.77928 

2.282677 

0.255888 

0.647773 

m 

PCOD  = 

436.7418 

2.07E-09 

3.683165 

419.3357 

N  = 

12.10356 

-8.4E-10 

14.80566 

5.867708 

Efil 

N02  = 

0.843553 

4.56E-09 

86.60573 

73.77928 

51 

Ox  Consu 

220.5524 

1037.525 

278.8705 

81.895 

5589.349 

52 

53 

E3 

SCOD  = 

0.677953 

-1.7E-10 

351.2799 

4.259734 

NO  = 

2.003096 

-9.6E-08 

106.7535 

0.749177 

0.483809 

1 .206775 

PCOD  = 

430.3321 

- 1 

8.06E-09 

6.861583 

565.8419 

l57 

N  = 

10.51041 

-2.4E-08 

117.8955 

12.44389 

8.035683 

m 

N02  = 

0.968254 

1.21E*07 

117.8955 

106.7535 

59 

Ox  Consu 

169.3293 

575.4189 

118.4963 

1973.154 

60 

61 

Ojr-4  /8 

' 

SCOD  = 

0.05163 

•5.4E-14 

55.17513 

0.786792 

-2.4E-09 

222.0987 

0.138376 

MEESmi 

mniiiii 

O 

PCOD  = 

417.8581 

1.21E-11 

15.66991 

1098.88 

Ea 

N  = 

7.637138 

-8.8E-10 

231.8302 

1 .958667 

39.00924 

16.07137 

N02  = 

1.077168 

3.27E-09 

231 .8302 

222.0987 

67 

Ox  Consu 

284.1862 

1 

808.703 

746.4932 

246.5295 

68 

69 

0x  5  /9 

70 

SCOD  = 

0.007291 

-2.5E-13 

3.904244 

0.057475 

71 

NO  = 

5.674778 

-4.5E-07 

111.9972 

0.010108 

0.587333 

1.401601 

72 

PCOD  = 

411.7101 

1 .96E-09 

7.969341 

541.3559 

73 

N  = 

6.230772 

•3.4E-07 

114.3294 

0.13866 

19.22638 

8.035683 

74 

N02  = 

1.10327 

7.88E-07 

114.3294 

1 1 1 .9972 

75 

Ox  Consu 

138.9334 

388.3726 

368.1407 

124.3169 

1173.107 

76 

77 

Ox-6/10 

78 

SCOD  = 

0.000373 

3.21E-15 

0.199943 

0.003082 

79 

NO  = 

'  9.230935 

-9.8E-08 

106.3925 

0.000542 

0.61499 

1 .403743 

80 

PCOD  = 

393.7959 

-1.6E-11 

7.981523 

517.8006 

81 

N  = 

2.39589 

-1. IE-07 

102.18 

0.007106 

18.40237 

8.035683 

82 

N02  = 

0.959743 

2.12E-07 

102.18 

1 06.3925  i 

83 

Ox  Consu 

129.6188 

374.661 

329.0196 

118.0957.  1114.459 

84 

85 

1 

86 

87 

_ ! _ J 

► 


82 


83 


84 


85 


|ee 

\fm 

Ibb 
Iqq 

BE 

|F77«1 

IF^ 

to 

to 
toi 
toi 
toi 


(concentration) 


(concentratio) 


1 

68.06068 

2 

65.65949 

3 

23.21862 

4 

21.82914 

5 

21.3406 

6 

4.657131 

7 

0.677953 

8 

0.05163 

9 

0.007291 

10 

1 7-Nov 


SCOD  Concentration 


BBl 

BSI 

mm 

BQI 

BE! 

BE! 


200  i\ 


S  150 
a> 

E 

u 

i  100 


BSI 

BE! 

BE! 

mm 

mm 

mm 

I^K. 


Calculated 


Plant  Data 
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NITRIFICATION 


Term  la) 


Term  |b) 


Term  (c) 


331  Term  (d) 


Stage 


Depletion  of  NH4  via  heterotrphic  uptake  of  soluble  COD 


Depletion  of  NH4  via  heterotrophic  uptake  of  particulate  COD 


Adding  NH4  by  death  of  heterotrophs 


Depletion  of  NH4  by  Nitrosomonas 


Term  lb) 


mg/L-day 


Term  la) 

mg/L-d 

mg/L-hr 

80.83992 

3.36833 

35.04241 

1 .4601 

iMHcv;«!;irancgri»j 

153.7224 

6.405101 

26.77256 

1.115523 

2.107 

0.087792 

0.298322 

0.01243 

0.015288 

0.000637 

mg/L-hr 


Term  (c) 


mg/L-d 


10.282 

0.428417 

3.739775 

0.155824 

2.197605 

0.091 567 

0.799622 

0.033318 

8.156236 

0.339843 

3.224434 

0.134351 

3.638512 

0.151605 

1 .438896 

0.059954 

1.293142 

0.053881 

0.511449 

0.02131 

43.62304 

1.817627 

17.28847 

0.720353 

43.12526 

1 .796886 

17.28847 

0.720353 

41.96347 

1 .748478 

1 7.28847 

0.720353 

41.36485 

1.723536 

17.28847 

0.720353 

17.28847 

0.720353 

255.173 


253.6478 


245.9755 


219.8365 


10.63221 


10.56866 


10.24898 


9.159852 
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Appendix  B:  Storage  Model 
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mi 


TABLE  A.1 


RT  = 


N  I  0 


(INPUT  VALUES  FOR  STORAGE  MODEL) 


5 


So  = 


Qras  = 


S(10)  = 


Qir(ax-1 1 


S9  = 


SI  = 


Y(h)  = 


Fdn  = 


Xh  = 


Ks  = 


Kos(het)  = 


01  = 


N01  = 


NOs  = 


VI  = 


NOo  = 


N010  = 


N09  = 


Kd(het)  = 


Xnb  = 
30|Kn(ns) 
20  A  = 


31 

0.5 

0.175874 


Fh  = 

0.1 

Kno  = 

0.75 

Um(H>  = 

2.8 

Kn(nb)  == 

0.6 

2977.459 1  Kon(nitrifi) 


20  Koh(inhib) 


0.1  Vltotal)  = 


K  (sto)  = 


0.75  Kss  (sto)  =1 
0.156|sr  = 

Ss(10)/Xh 


10iSs(9)/Xh 
Xltotal)  = 

0.1758741  I 


0.06 


0.937 


0.002 

30 

0.000115 


0.000388 

3910 


Qir(ax-3)  = 

40 

06  = 

3 

V2  = 

0.145 

V3  = 

0.5047 

V4  = 

0.124 

V5  = 

0.124 

V6  = 

0.3394 

V7  = 

0.4648 

V8  = 

0.9296 

V9  = 

0.4648 

V10  = 

0.4648 

solve  for  X(h)  and  P(10) 


X(h)  and  the  initially  assumed 


P(10)  value  are  estimated  by 


X(h)  =  2977.459 


PI10)  =  392.4781 


clear  all  concentrations 


4.55E-13  411.3998 
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C  I  D 


i  I  PHOENIX  PROJECT!  \  (STORAi 

MoiM  FOB  ANAL  YSiS  (based  on  ST0RQE13.XLS} 

3  1  Q(o/  -  20;Q(fas) = 9. 3S;Xi(ol = 20  Xnsand  Xnb  are  fixed 


AX'! 


6  ISCOD  =  I  54.07295 


NO  =  2.010418 


8  IPCOD  =  479.3755 


9  N  =  21.20262 


N02  =  0.537893 


1 1  lOx  co/jsu 


12|Ss/Xh  = 


13 


14  ISCOD  = 


15  NO 


16  IPCOD  = 


F  I  G 


(STORAGE  MODEL} 


SCOD  = 

30 

NO  = 

2.164613 

PCOD  = 

439.3013 

N  = 

14.23665 

N02  = 


27  I  Ox  Consu 


28  Ss/Xh  = 


SCOD  = 


31  |NO  = 


PCOD  = 


N  = 


0.012494 


AX-A 


30 


1.804206 


438.6728 


14.18234 


PCOD  = 


41  IN  = 


N02  - 


Ox  Consu 


44 


SCOD  = 


47  |N0  = 


PCOD  = 


QX:ZJ& 


30 


2.261081 


433.3695 


7.05E-09 


4.9943351 

2.786222 

1.0148531 

:  0.307213 

-1.7E-07  79.2693  16.62348  42.32849 


-3E-07  476.8673 


0.299933 


15.77506! 

8.423622! 

3.34811 

!  0.299933 

1.31E-07 

1  18.44067 

3.884555 

2.56E-09 

1  o; 

3.669804 

1  1.965605 

0.782382! 

1  0.283993 

73.831041 

1.647538 

5.207603 

416.0978 

0.91 5883  i  0.00927 
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miu. 


02  = 


2.89E-13 

1 

-3.6E-13 

51 

Ox  Consu 

103. 1783 

52 

Ss/Xh  = 

0.008919 

NO  = 


30 


3.411185 


426.9998 


10.95301 


0.970871 


Ox  Consu 

i  141.2656 

Ss/Xh  = 

\  0.005466 

Ox-4  78 

SCOD  = 

30 

NO  = 

5.822783 

N  = 

7.687354 

N02  = 

1.079174 

Ss/Xh  = 


SCOD  = 


NO  = 


PCOD  = 


N  = 


76  ISs/Xh  = 


77 


SCOD  = 


NO  = 


PCOD  = 


0.00147 


31 


3 


7.046403 


408.5141 


6.184886 


0.00051 


2  = 


Ox  Consu 


Ss/Xh  = 


30 


10.5886 


390.7364 


2.32638 


0.952856 


128.9668 


6.42E-05 


G 


14.74569 


86.784851  73.831041 
279.44721  81.952461  1226.721 


5.867708  0.740685 


-3.3E-08I  106.86371  2.2076261 
~6.2E-08|  7.67107V  561.4603! 


-6E-14  118.1122  32.155461 


-1.4E-13  118.1122  106.8637 


380.3214  118.6187 


0.5451071  1.349143  0.00907 


19.9196]  8.035683;  0.663857 


1208.625 


-4. IE-081  222.24641  2.760602|  1.1781621  2.831391 1  0.005671 


-9.7E-08 


1.6E-121  231.92331  37.2083  38.72588  16.0713710.384087 


222.2464 


776.8415\  746.7929  246.6935  1173.313 


0.60099 


-3.6E-09 


-1. IE-08  8.091338  537.1535 


6.75E-14  114.2664  8.9311781  19.083571  8.03568310.184386 


-1.8E-13  114.2664  112.0305 


384.338  \  367.938  124.3538 


-1.2E-07  106.09751  0.10677  0.6209731  1.40638810.000439 


-4.5E-07  7.996558  513.7777 


18.2621  8.03568310.028184 


101.6556 


372.5435  \  327.3309  117.7682  1105.402  ! 


-1.6E-07  0.0131 


94 


95 


96 


97 


98 


99 


ES 


EB 

Ieb 

\m 

|eb 

|eb 

|eb 

|g!«ni 

|eq 

Ieb 

\m 


ODI 

QQI 

EIEI 


mm 

EiBI 


mm 

eebb 


Calculated  Oxygen 


E]^l 

Fmm 

mm 

iMm 

mm 

EIBI 


E£B| 

m 

eb 

EB 

m 
m 

E3 
122 
E2D 
ES 


141.2656 


279.2338 


138.2708 


128.9668 


Calculated  Oxygen 


80.30819 


80.28857 


79.35155 


78.58649 


73.29851 


E 

I  40 
o  30 


100 


102 


cz 

A 

B 

C 

_ D _ 

E 

P 

G 

H  1 

□E 

li'IW  H  III! 

CE2 

Term  (a) 

Depletion  of  NH4  by  heterotrophic  utilization  of  soluble  degradable  COD  in  cell  storage 

Term  |b) 

Depletion  of  NH4  via  heterotrophic  uptake  of  particulate  COD 

^  -  ■ 

: 

EES 

Term  {c) 

Adding  NH4  by  death  of  heterotrophs 

-  ,  . . . 

CS] 

Term  (d) 

Depletion  of  NH4  by  Nitrosomonas 

- ^  . . 

gygi 

Staae 

Term  (a) 

Term  (bl 

Term  (c) 

Term  td) 

mg/L-d 

mg/L-hr 

mg/L-day 

mg/L-hr 

mg/L-d 

mg/L-hr 

mg/L-d 

mg/L-hr 

FFn 

1 

32.01497 

1.333957 

17.8604 

0.744183 

6.505469 

0.271061 

0 

0 

gFW 

2 

27.79585 

1.15816 

15.46255 

0.644273 

5.647369 

0.235307 

0 

0 

EED 

3 

16.69036 

0.695431 

6.633862 

0.276411 

0 

0 

EE 

4 

15.85165 

0.660486 

6.309534 

0.262897 

0 

0 

5 

14.83133 

0.617972 

5.911334 

0.246306 

0 

0 

EEi] 

6 

77.18755 

3.216148 

43.44634 

1.810264 

17.28847 

0.720353 

255.7008 

10.6542 

EE 

7 

69.18127 

2.882553 

42.85627 

1,785678 

17.28847 

0.720353 

EE 

8 

40.02614 

1.667756 

41.65865 

1.735777 

17.28847 

0.720353 

249.4872 

10.3953 

EB 

9 

19.2151 

0.800629 

41.05759 

1.710733 

17.28847 

0.720353 

gkn 

10 

2.937111 

0.12238 

39.29023 

1 .637093 

17.28847 

0.720353 

218.7082 

9.112841 

EIB 

. 

E& 

EE 

ee 

Total  Nitrification 

ee 

441 

Staae 

ma/L-d 

ma/L-hr 

- 1 

EE 

1 

56.38083 

2.349201 

443 

2 

48.90576 

2.03774 

444 

3 

54.58<,o4 

2.274189 

445 

4 

51.75638 

2.156516 

- 1 

446 

5 

48.3474 

2.014475 

447 

6 

393.6232 

16.40096 

448 

7 

383.4401 

15.97667 

449 

8 

348.4604 

14.51918 

450 

9 

323.4012 

13.47505 

451 

10 

278.224 

11.59267 

Appendix  C:  List  of  Mass  Balance  Equations 
for  Storage  Model  and 
Michaelis-Menton  Model 
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Mass  Balance  Equations  for 
Michaelis-Menton  Mode! 


Ammonia  (NHLd-N)  Mass  Balances 


Stage  1 


Mm(NS} 

(  o,  \ 

Y 

V  ^NS  y 

(K,+N,) 

-n,[(a.s'.)+(fl„.5.,)+(a<.,s,)-(a+e«.+a<,;X^.)](^-,) 


-r. 


H 


(Q.P,)HQc  +  QfUS 


VA. 
e. 


+(0(1, 


»  j 


-4[-(a+eK.5+a(i)K^)](/v) 


HK,)iX„)(VrKF,)(Fa.) 


(  NO, 


Stage! 

i"^=o=(a+e«-,  +a<.,)(Af.)-(a+a«.  *Qh«xk) 


(  Mm(.NS) 

(^.sX/^2) 

(  0,  ] 

l  j 

(K,+N,) 

[Kos+oJ 

[(Go  QfUS  Qtr{\) K^l )  ~  iQo  ■*■  Qtus  GmD X^l  )](^W  ) 


105 


-^[(a+ews  +a,„)(/i)-[(a+e»,,  +a<„x^,)]](Jv) 


+(«^KJf,)(n)<^»X/w» 


NO^ 


VM?,  +  iVO, 


Stage  3 

5 = ^ = + fi„.)(^,) + (a,.)(  A'.)  -  (a + G„  +  a„ + a„x  a^.) 


/^iw(yj)  X^3 )  Qa 

Yss  j{K,^N,){K,,^oJ^ 


) 


^h[(2o‘^G|MS  (2o'^Q*45  "^OhI) '^Qlr<J)K'^3)](^»r) 

-n,[(a+aus  +a<.,x^5)+(a,„KP.)-[(a+au.  +a,.>+a,»x^’,)]](^',) 


4<A',KJr„x>',X^vXf, 


ON 


i  NO,  ^ 
\no,-^no,} 


Stage  4 

^ = 0 = (a + +  a,,, + a<.,KJv. )  -  (a + +  a,„ + ao,x^.) 


Mm(.ffS) 

f  Q.  \ 

<  Yns  > 

iK,+N,) 

^ff[(fio  '^0IMS  '^Qir{l)  "^QkaS  "^Otril)  fir(3) K54)](^^ ) 

“^ff[(!3o  "^QlUS  "^QiriX)  ■*’i2tr(3)K'^)“[(Co  "^QraS  "^QiriX)  ■^CL<3))(-^)]j(^  ) 

S  NO,  \ 


HK,)(X„W.WsWvs^ 


NO.-k-NO, 


106 


Stage  5 

dN 

y-^=o=(a+e^s  +aKi)+ao,)(^4)-(a>+ei«.  +cu,>+a<.,)(^.) 

_  o,  . 

1,  y,,  )(k^+n,)[k,,+oJ 

[(fio  QmS  ■*■  Glrfl)  ■*■  Qlrii)  )(^4  )~iQo'^  QkaS  ■*■  Qr(l)  ■*"  Qhii)  )(^i  )]('^y  ) 
~^ff\(Qa  QkAS  "^GUl) '^QiK3)K-^)“[(Q)  ■*'SpUJ  ■*'Qr(l)  ■*'flr(J)X^)]](-^) 

\  y 

Stage  6 

dN 

^  = 0 = (a + +  a^i, + a<3,)(^5)  -  (a + e*..  +  cui, + q^.,)w 

I,  J(/:^+iv.)Uoix+a/ 

~^[(a  ■^aM5  ■^ar(3)K^5)”(a  ■^■GjMJ  ■^'SkI)  ■*'Sr(3)K‘^6)](-^  ) 

~^[(a'*‘2iM5  ■*'Q(r<l)  ■*‘a«<3)K'^)“'(a  ■'■fiMS  ■*'2^(1)  ■^ar(3)K^)](-^Af  ) 


0  =  (Qi+QluS  +Qr(l)+Q,rii))(^6)'~(Qe’^QlUS  +  Qr(l)  +Qr<3))(^7) 

_  f —2i—\v ) 

"^QraS  "^Oiril)  2lr<j)  "^QrAS  '*‘Qr(J)  Qr(J)  )](^^  ) 

[(fio  fltJS  "^Qk!) ■*'Qr(l)  ■^i^3)K'^)](^  ) 

0  =  (Co +0145  +Qr(l)‘*'CU3))(^7)-(Co+CjM5  +  C^J)  ■•■  CUs) X ^» ) 

_  O.  \j. . 

~^[(fi)  "^QlUS  ■^QkI)  ■*'0r(3)K^7)~(Q>  "^ClMS  OnD  ■*■  CK3)K'^i)]('^  ) 

~^[(Q)  ■*"ClMS  "^Oird)  ■^0r<3)K-^)~(C<)  QkAS  Qr<I>  ■*"  CL(3)X^)](^  ) 

■K/f««)(-!f,)(»'.)(^',)fT^l 


Stage  9 


= 0 = (Q, + +  Oo) + -(Qo+ Q^s  +  ao) + ao,)(^,) 

_ (^«)W f  o.  . 

I,  Y,,  J(a:,+a^,)Uo^+^J 

+Q,<1)'*’2U3))(‘^|)~(2o'*'2|UJ  +Qr(»)  ■•■QK3))('S'9)](^Ar) 

~^[(0O  "^QllAS  "^Oirm  ■^flr(J))(^)~(2o  '^'QrIS  i^3)K^)]('^  ) 

■K3r«.,)(Jr,)(K,)(/-47-^ 

Stage  10 

f'^=o=(a+e«,)W)-(a+a„)w.) 

r ('^^wyK^m)  Q|n  y 

I  i;.  J(/f„+w„)U„+o../  "•' 

-j'.Ka +a«y)(s,)-(ft +fi„,)(s..)Kf,) 

-j;[(a +e«.)(/’,)-(a  +e«.)(^.)H) 

-KJf«.,)(3r„)(K„)(/v)|7^ 

\*«r+^ioy 


109 


Stage  1 

>'^=o=(a)(w.)+(0u.)(oi)(^o,.)+(a,„)(Aro,)-(a+a«+a,.,)(w,) 


110 


Stages 

® "  (a +e_+a,.,)(  w,) +{a,„)(wo,)  -(a +a„,  +a,„)(Jvo.) 


i-i.42y„ 


2.86 


V  4  y 


zw 


V 


jA 


NO.  +  NO, 


iJ 


1.42 

2.86 


Y  m 

V  ^<1*  ■*"  ^3  >  V.  ■*■  ■^*^3  y 


"a* 


v^M+^3yv 


m 


NO.  +  NO, 


iJ 


Stage  4 

y^=o=(Q, +s.+a<„  +a<.,XAro,)-(a  +&.,  +a,„ +a„,)(w.) 


ri-i.42r,Y//„„l 


2.86  A  n, 

1.42 


onr 


v^,+‘S'4yv 


.  ^4*  ■*■  ^4 


jUo,  +  A^oJ 


2.86 


f  k  ^ 

f  1 

\NO,  +  NO^^ 

^i-i.42y„ 
,  2.86 


NO^ 


k^+Oj(NO^  +  NO, 


K 


stages 


dNO, 


= 0 = (a + + 0(1) + a^3,)(  -  (0, + Q^s  +  Ok,) + Qk3,)(  ^^3) 


2.86  )[  Y„  j 

1.42 


Mb 


V 


K.  +  S. 


NO. 


sA"-** 


k^-^O.  [NO.  +  NO. 


2.86 


K  ] 

f  1 

<  ^«ik  ■*■  ^5  J 

[iVO,  +  iV(2j 

1-1.42/, 


2.86 


^KiKM)xJ 


f  *-  1 

(  1 

\k^+0,^ 

[  NO.  +  NO,) 

n 


Stage  6 

^^=o=(a+Q^,  +0(3) +i2U3))W)-(a +0^1) +0^3, )(^^,) 


f iw®) Y (^^)(ArQ2,)Y  0,  Y 

i;*  A  k,^n02,  Xk„+oJ* 


fjS^Y-^Y— 

I  2.86  A  Yb  j  °"[x,  +  S,Xk^^O,XNO^-^NOj 


L42 

2.86 


(KXX^XFa,) 


NO. 


2.86 


^oh 


NO. 


NO.  NO. 
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Stage  7 

y^=o-(Q,+Q^s  +Q»»+&i„)(No,)-{Q,+Q^,  +a,„+ao,)(w,) 


+ 

V 


MmiNB) 


Y 


(X,,)(N02,) 
+  N02, 


JU-.+OJ 


n-i.42y„ 


V 


2.86 


V  4  y 


m 


K^+S. 


V 


\r 


7  A 


k^+a 


NO, 


’7  A 


A^o +m 


\v., 


1.42 

2.86 


V 


A 


m  +  NO. 


Tj 


ri-i.42y„^ 

2.86 


^7 


ATQ, 


yv 


NO,^NO,j 


% 


Stages 

f'^=o= (a  +e«..  +a«- +a,„)(/^o,)-(a +e«..  +a,„ +a<.,)(^a) 


(Jf™)(JV02.)'| 

(  o,  ^ 

1  y.B  J 

t  K,+N02,  J 

1-1.427^ 

2.86 


\( 


A  y. ) 


DN 


s,x„ 


yv 


.  ^4* 


ATO, 


yv 


NO^  +  NO, 


1.42 

2.86 


iK)iX„){F^) 


\r 


NO, 


{k^+0,){NO,^NO,) 
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i-i.42r„ 

2.86 


\ 

(  k  Y 

y^iK,)(p,)x„ 

^oh  . 

<k.^+oJ^ 

NO, 

NO,  +  NO,, 


Stage  9 


^  j  ® = (a + e,„s + a,,, + OoiX^o.)  -(&+ +  a,., + a<,,)( m) 

at 


NmjNB) 


V 


V  A 


(X^)(N02,) 

K^+N02^ 


V 


a 


W 


A-i.42r„^ 

Fm(H) 

f 

{ s,x,  ] 

1  2.86  ) 

[  y»  j 

-*  DN 

[k+sJ 

m  \ 


1.42 

2.86 


(k,)(X„)(F^) 


(  k  \ 

"'ah 

[no,+noJ 

V 


i-i-42y„ 

2.86 


F„iK,)(P,)X„ 


r 


V 


*ch 


NO. 


A 


NO^  +  NO, 


K 


Stage  10 


dt 


0 = (a + Qn.s  )iNo,)-  (a )( Aro„) 


'hb  a 


K,+N02,, 


A 


a 


10 


ife_+a 


10 


ri-i.42y„^ 

*  DN 

S,A  I 

(  *-  1 

r  NO,,  A 

V  2.86  ) 

1  4  J 

1  ^o*  **■  Ao  / 

[no.+no,J 

1.42 

2.86 


(k,){X„)iF^) 


V 


kph  +  ^10 


A 


NO, 


10 


NO.  +  NO, 


10 


10  J 
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fl-1424 

I  2.86 


f 


_^Y 


V 

no,-¥NO,,} 


Nitrite  fN07'*l  Mass  Balances 


Stage  1 

= 0 =(axAf02.)+(a„,  KAf02„)+(a„,xAfO2,)-(a +a«,  +a<„)(Af02,) 


Stage! 

= 0 = (a + a_ + a<„)(A'02,)  -  (a + q„,  +  a,.,xJvo2,) 


Stage  3 

f'^^=«=(a+fi.+a,.,KW2.)+(a,„XJ'02,)-(a+a«+a,.,+a,»XAf02,) 


Stage  4 


ydN02, 

dt 


0  =  (Si.  +  Qm,  +  Qkd  ■*■  Qr<3))(^^23)  -  (Q,  +  032^5  +  0^,)  +  0^j))(NO24) 


stages 

dt 


= 0 = (0 + o. + + aK„)(^^o2,)  -  (0 + 0,^,  +  0^, + 0^3))(^o2,) 


115 


Stage  6 


= 0 = (a + + a,,, + a<„)(  wo2,)-(a +at,. + + a„,)(A'02.) 

Jf^asiX  ■yy.Af.  Y  o,  Y,, 
t  JU„+a'.Aa:.+o.J'  ‘ 


Mm(NB)  Y 


^NB  Jv  ^N(m) 


N02.  Y  Q._  V> 


Stage  7 

= 0 = (a + e«.. + a,., + a,»)(wo2.)  -  (a + e«., + a™ + Oo.ka'o^,) 


if  ^ss^l  Y  ^7 


Mm^rn)  Y 


^SB  Jl  ^N(ffB)  ^N02,Xk^+0, 


Stages 

= 0 = (a + + a<„ + Qh«)(ho2,)  -  (a + e*„ + a«„ + q^^mnoi.) 

at 


+f^T#^T^7%rV.) 


y„  Xk,+n,Xk,+o, 


Ysb  JU^(™)  +  ^^2.A^^+o, 
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Stage  9 

V  = 0 = (a + + a,,, + a,,, )( JV02, )  -  (Q, + + a,,, + a,., )( ^^02, ) 


^m(NS)  Y  ^NS^9 


yss 


y.B  Ak,,,,,  +  N02,j[_k^-*-0, 


Stage  10 


= 0 = (a +e»,»  )(wo2,)  -(a +fi„,  xnoij 


^m(NS) 


Y  ¥  0,. 

A^iV  A^«.+^k 


X^„N02 


A  ^v(Mr)  A  +  O, 


SCOP  Mass  Balances 


Stage  1 

=0 = as.  +e*.,s.. +a,„s,  -(a  +e««  +a,„)(s, )-(«.)(»",) 

at 


where;  /l,  =  A^(^av) 


^,+5,A^^+o,A^o.+Aro. 
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Stage  2 

= 0 = (a +0. +a„.)(s,)  -(a +e,u.  +a^„)(s,)  -(«,)(n) 


where: 


f  c  y 

=  K(F^) 


+  A 


m 


m+m 


2/ 


Stage  3 

K^=o=(a+a-+a,.,)(5.)+(a<,,)(s,)-(a+au,+a,.,)(5.)-(^)(f',) 


/ 


where:  =  A:(F„)| 


/ 


m  m 


3j 


Stage  4 

»"^=o=(a+e„+a,„ +a,.>K'^.)-(a +e«. +a,„ +a,.,)('5.)-(0(n) 


where:  =  K{F^)\ 


s,x„ 


V 


'ah 


NO, 


4  J 


NO,  +  NO, 


Stage  5 

»'^=o=(a+a.+a,.,+a„,)(s.)-(a+a.*+a,„+a„,)(5,)-(^.Kn) 


where:  R„=^KiF^) 


Mil 


\r 


A 


m 


NO,  -I-  A/O, 


Stage  6 
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Stage  9 

>'^=o=(a+e«.,+a,„+a<«)w-(a+e...+a,.,+a„,)(5.) 


rA:s,x„'| 

f  1 

U,  +  5j 

{0,+k„} 

r. 


(ks,x„^'^ 


k+S. 


V 


9  A 


m 


A 


NO,  +  NO. 


K 


9  y 


Stage  10 

^ = 0 = (a + )(s, )  -  (a + e,^  xs,. ) 


*^10 


a 


10 


yv 


^10+^. 


KS^ 


ifc +5. 


10  A'^d, 


k.^+0. 


NO. 


iS- 


'wyv 


NO,  +  NO,,j 


10 


PCOD  Mass  Balances 


Stage  1 


f^^=o=(ax/»,)+(a+0««K/i.)-^^^^+(ai.,xp.)-(a+0««+a<,,x/;)-/>r,) 


where:  r=F^(KPX„)\-^^^^ 


k^+Oj[NO,  +  NO,^ 


Stage  2 


f'^=o=(a+a.+a<.,)(/5)-(a+fi.u.+a,„)w-'-,(f;) 


where:  =  F„{K^P,X^) 


V  ^2 


NO, 


NO.  +  NO. 


Stage  3 

>"^=o=(a+Q.+a,., )(^.)+(e«,x/’, )-(&+&« +2-0, 


where:  r^  =  F^{K^P,X„) 


\r 


''eh 


A 


NO. 


NO+NO 


ij 


Stage  4 

o=(a +s.+a,„ +a„,)(/’,)-(a +e,t„  +a,., +a,„x/’.)-'-,(n) 


where:  =  F^{K^P^X„) 


k^+O^ 


A 


NO^ 


NO,  +  NO^ 


Stages 

f'^=o=(a+e«.+a,.,+a«.,x/’4)-(a+e,„.+a<„+a<.,x/;)-';(f;) 


where:  r  =  F„(*r 


V' 


NO. 


k^+Oj{NO,  +  NO, 
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Stage  6 

0= (a  +Q^  +a,.,  +aK.,x^.)-(a  +eu  +a,., +a<.,)(/’.) 


<krx„v,) 


k^  +  O. 


w 


Stage  7 


= 0 = (a +0^.) +a.3,)(^*)-(a,  +(?»..  +a<i, +(2k3))(^7) 


A 


NO, 


NO,  +  NO, 


<K^P,X„V,) 


^  O,  ^ 


Stages 

^  ^  ~  "^Qius  '^Qiril)  ■*'Gh^J)X^)”(2o  '^QkaS  "*■  6lf<3) X ) 

at 


-{KRX„V,) 


\f 


+0, 


NO, 


\A 


NO^  +  NO, 


122 


-(KAXM 


'  O.  ^ 

K+o, 


Stage  9 

dP 

= 0 = (0, +0(1) +a<3))W'(a +a<.) +a(3,)(n) 


-{krx„vA 


NO^\ 


~{KRX„V,) 


^  O.  ^ 

V  ^a#  +  ^9  y 


Stage  10 

f'^=o=(a+e«,)(^.)-(a+s,„,)W 


\r 


<KR,X„vA 


^th  ■*■  ^10 


m 


10 


NO.  +  NO, 


10  / 


^K,P„X,vA-^  ' 

+C/,o  y 


123 


Oxygen  Coiismw|itioii 


Stage  6 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  [(^)(r,  +  1.42(i:,)(Jf„)](^'.)+(3,22)(v,„,)(K,)  +  (l.n)(r«,,,„,)(r.) 

wha%; 

^  (ks,x„Y  o,  \ 

.4  =  l-1.42(y„) 

^ns  KK,^N,Xo,+k^) 


(  0,  \ 

•nOjin)  -  Y 

\  J 

{  K,+N02,  ) 

Stage? 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  [(>«)(n.  +  1.42(/:J(^J](F,) +(3.22)(r^„^,„,)(K,) +(1J  l)(r^^^,„,)(K,) 


^ere; 


r_  = 


KS,X„ 

k^S, 


O. 


A 


A+A:. 
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^  =  1-L42(4) 


’NHtitt*} 


‘y»  / 


O, 


’NO^iax) 


I  y.n  J 


(X,,)(N02,)Y  O, 
K,+N02,  + 


Stages 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 


=  [(^)('',  +  1.42(i:,)(X,)](K,)+(3.22)(r^.,„, 


KK.)+(l.lIXr^,„){K,) 


where; 

^  _  (ks,x„  Y  o,  ] 


.4  =  1-L42(y„) 


MmjNS) 

V  J\ 


xm)Y  o. 


Mm(m) 


(X^XNOl,) 


A  k,+no2. 


O, 


o,^K 


Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  +  1.42(*:,)(  jr„)](K,) +(3.22Kr™.,„)(K,)  +  (tl !)(»; 


^^ilere: 


K^]  O, 


^,+*^9  A^9 


^  =  1-1.42(4) 


126 


Stage  10 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 


^^lilere: 


/■_  = 


K  *^10 


a 


10 


A 


^  =  1-L42(y„) 


\  ^ra  J 


x,ANJ 

\^N  J 


o. 


10 


V^io  "^^0 


V 


(X^)iN02J]  Q 


\  *1*8  A  *^N 


K„+N02 


10 


10 


127 


Storage  Model  Mass  Balance  Equations 


(Mass  balance  equations  for  PCOD  and  Nitrite  in  the  Storage  Model 
are  identical  to  the  mass  balance  equations  listed  previously  in  the 
Michaelis-Menton  Model) 


Sc/Xnr  Mass  Balances 


Stage  1 


+f^W,  -ft5.)(0.7)-(/e.^)(K,) 


vdiere; 


^5,V  ~  ^jfo 


(  *0*  ] 

(  NO,  ] 

[no,+noJ 

Stage  2 


y--  = (a +0«+aKofe«)/^/r)-(eo 


^  1  ^ 


where; 


128 


S,V  M>  f  ,  O  /  ir  I.  I  /->  ■  X 


k„  +  A  +  O,  jl,  NO,  +  NO, 


(/V) 


Stage  3 


^  . . -  (a  +Qk3)K^^3)/-^J 


+('^l[(a+fi,+a,,,)W)-(a+fi.+a,,,)(s.)](i.<>)-K„)(f'.) 

\^H  J 


v^ere: 


/?  =Jt  ^S(3)I^H  ^^3  fp  \ 

"U.+WJ^JU-+oJUo,+^'o,J  “' 


Stage  4 


:i 

^=(a+a.+a.„+< 


^  a  =  (a + a« + a.., + ax3,)(5.,3,/^„)  -  (a + Sxo + +  aK3,)(‘5.(.)/^K) 


-M(K) 


v^ere: 


R  =ik  ^S(4)I^H  ^ck  ^Qt  fp  \ 

’‘\K^S„,,lX„Xk,,^0,lN0,^N0j 


129 


Stages 


= (a + s. + a,,, + Qjis,jx„)-(Q, + a  Kl)  Qras 


where; 


=  k. 


\r 


\r 


'■i* 


K  +  S^JX„  jl,  +  O,  A  NO,  H-  NO,  J 


NO, 


(Aw) 


Stage  6 


V 


dt 


~  {Qo  ■*■  Qnu  **■  Qirin  (Qi  Q/(l)  ■*“  Ql*S  ■*"  SiK3))(A(6)/'^h) 


-K- )(»">) 


where: 


A(6)/^g  r  0}  ^ 

A»  ^  V^at  ■*■^6  / 


Stage  7 


^1  -S(7) 


= (a +a„+aK« +Qj(s,jx„)-(<i+a„.,^Q^ +q^Js,jx„) 


-(kAk) 


wiiere; 


As,:/  ~  ^ 


''  S,„,/X„  Y  o,  ^ 

As(7)  /'^H  ;  I,,  ^0.  ^7  y 


130 


Stages 


^  ^  ^  — = (a +a.3,)ten/jrj-(a +aJ{Ssj^.) 


where: 


^s,v  =  * 


O, 


j\ 


Stage  9 

=  (a + a. + s.„ + Q„)te™/-y„)  -(a + e«, + +  s,J(s«/-^„) 


w^ere: 


*^5(9)  ! 


O, 


Stage  10 


V  ^  X— ^ = (a +eJ(s„,/^,)-(a +e«-J(s,„.,/jf„) 


^S.I/  ~  ^ 


a 


10 


^^1le^e: 
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SCOP  Mass  Balances  (StoraEe  Model 


Stage  1 


H  0= as.  +{&,«)is,)HQ^)(S„)-(&+&,„  +Su,)(S,) 


-(as.-as,)(o.70) 


Stage  2 


n j= 0= (a +e«, +G,„)(s,)-(a +0,1, +a.»)(s;) 


-[(a +a.+a,.,)(s.)-{a +a.+a,„)(s.)](on)) 


Stages 


''l^-^j=o=(a+e«»+a<„xs.)+(a<„)(s,)-(a+a,.,+au«+a«,)(s,) 

-[(a + a. + a<.,)(s, )  -  (a + 2- + a<.,  xs.  )](i.o) 


S3  =  S4  =  S5  *  Sfi  =  §7  =  Sg  =  S9  =  Sio 
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Ammitnia  fNHj.Nl  Mass  Balances 

Stage  1 

f'^=o=(axjv.)+(e^,)(w„)+(a,„Kw,)-(a+e,„,+a<„)(w,) 
(,  r„  j(K,+N,}[K„+oJ  ’’ 

(v P  \ 

-Yh  (a/*o)+(a+0iu.)(^o)-I^J+(a<„x/»,)  (/v) 

-Y„[-{Qo+Q^s'^Qn.)W)\Fs) 


HK,){X„W.W,W^)\ 


] 

NO,+NOj 


wiiere: 


Stage! 

'"^=o=(a+e«„  +a<„)(A',)-(a+a..  +q»„)w 

f  o.  \r ) 

I,  r„  j(K,+Nj[K„*0j  ’’ 

-y«WFs^WMK) 

-j'.Ica+Su.  +a,.,)(^)-[(a+au.  +a.,x^,)]](f») 


133 


where; 


D  ^afc  /p  \ 

^\k„HS^JX,)j[k^^0j[N0,  +  N0j'^  ^ 


Stage  3 


y^=o=(a+Q^+ a..,){N,)Ha,j{N,)-(Q, + g„ + g,, + g„)(  a^,) 


Mm(NS)  {X, 

Yhs  )(<K, 


NS  X^a)  1/f/ 1 


-4[(q,+0r.5  +Q^„K/*2)+(a,3,xp,)-[(a+eiu.  +a<.)+a(3)X^3)]](/v) 


A^O,  +  M?, 


where: 


R  =k  ^s(3)/^H  I  k^  I  [/ p-  \ 
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Stage  4 

dN 

= 0 = (a + +  Oki, + a(3))(^,)  -  (0, + Q^s  +  aa, + 


f^miSS) 

iX,sXN,) 

i  O4  \ 

<  YnS  j 

iK,+NA 

-};[(a+0«..  +a(.)+a(i)K^3)-[(a+e/u,  +a<„+a<3,K/i)]](/v) 


+(/:,)( ;r^)(f;K/vKiw)|^ 


NO,  ^ 
ArO,  +  M?J 


where; 


Nsju  ~  ^ao 


SsJX„ 


(  k. 


\r 


k„  +  (5^4, /^«)  A  A  ^O,  +  ATO,  j 


NO^ 


(fa.) 


Stages 

dN 

*"^=o=(a+a«.  +a<„+ao,)W)-(a+e„,  +a,„+a„,)(Af,) 

at 


f  Mm(.5) 

iX,s)(N,) 

(  0,  ^ 

1  ^AfS  > 

{K,+N,) 

m) 


-YA^H(fs^)(mf.) 
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( 


NO.  \ 
M3,  +  ArC>J 


v^ere: 


R 


•SM 


V 


K  +  ){  k^+Oj[  NO,  +  NO, 


NO, 


(Fa.) 


Stage  6 

dN 

=o=(a+i2^  +a,„+a,»)W)-(a+au.  +a,„+a<.,)w 


(^iV^X^o) 

(  o,  \ 

Y 

\  *NS 

(N,+N,) 

[Kss+oJ 

-»'»((a+  '^0^3)X'^)  (Go  ■*■0115  ^CUl) ) 


k_+o. 


'•sju 


=  K 


K 


^S(6)I^H  ^  Os  ^ 

+(‘y^«/^«)JU«+oJ 


where: 
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Stage? 

= 0 = (0, + Ou.  +  a,,, + cu.iKA',)  -  (a +Q^  +  a,., + a,„)(Af,) 

(  j(A:,+Jv,)U»+o,/ 

“^H[(fio  "^QraS  "^OiriX)  Qr(J)  K  A  )](^V  ) 

-KK^»XX„)(K)(F,{j^] 


where: 


^5,V  ~ 


SsaJ^H  Y  O. 


Stages 

^ = 0 = (0 + Ou,  +  Gki, + Q^«)(N,)  -  (0 + 0u,  +  0,1, +a,.,xjv.) 

i  ^■<».i'|(-y...)(JV.)f  _3l_  V ) 

I  Y,,  )(.k.*n,)[k.„*oJ 

-YA^M^xnmK) 

~^/f[(^  "^QlUS  "^QiriX)  ■^fir(3))(^)~(fio  "^QkaS  ‘^fir(I)  fir< J) ^ ) 
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where: 


^SJU  ~  ^jt» 


\r 


o. 


Stage  9 

dN 

v-^ = 0 = (G. +e«.,  +a<„ +a<s,)(^.)-(a +e^,  +aa, +a<3,)(^,) 


fu  ^ 

^m(NS) 

(^v,)(^,) 

(  o,  ^ 

<  Yns  y 

(K,+N,) 

(n) 


-YAXARs^)iK)Ws) 


"^QraS  ’^QhD  ■*'fiK3))(^)]('^) 


O, 


k^+a 


R-SJU  ~  ^mo 


V 


O, 


wiiere: 
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Stage  10 

= 0 = (a +at..  )(Af,) -ca +0t„  XA'„) 


NmiSS) 

f  0  \ 

'^10 

<  ^ffS  J 

<  ^10  > 

M^ere: 


^S,t/  ”  ^tu 


■*■^10  y 

f  ^  i  y  V 

^  +  (^j(IO)/'^w)  ;  ^10  J 


a 


10 


Nitrate  (NOi-N)  Mass  Balances 


Stage  1 


°  '  (S)  w) + (eu  )(«•  i)(A^o») + (a„,)(  w.)  -  (a + a«  +  aoiKw,) 


L42 

2.86 


/ — -  -  Y  — V, 

'  "  ""\k^-^0,XNO,  +  NO,P 


1-1.42J;^ 


2.86  ; 


F^(K,)iP,)X„ 


V 


NO, 


NO^  +  NO, 
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where: 


NO^ 

NO,  +  iVOj 
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Stage  3 

»"^=o=(ft+e>.+a„,)(Afo=)+(a„,)(A'o.)-(a+G,„.+a<„)(A'o.) 


1.42 

2.86 


f  1 

f  1 

(NO,  +  NOj 

1-1.427. 


2.86 


*-  1 

f  1 

Aw. +  <^3; 

wiiere: 


^5,U  *“  ^4»o 


^50)1 


k^HSsJJ^n) 


^  -  Wos) 


NO,  +  NO 


j/ 


Stage  4 

>"^=o=(a+G.+a,„+a,.,)(w.)-(a+a„,+a,.,+a<.,)W) 


1.42 

2.86 


(A:,)(^^)(F^) 


m 


NO,  +  NO^ 


1-1.427 


2.86 


*«* 


K 


M7, 


'^4 
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where: 


^ao 


^S(4)/-^H  Y  NO^ 

K  HSsjx„)  JU^ +0.  km + m 


iFv.) 


Stage  S 


dNO. 


>'~^o=(& +a,+a,„ +a«,)(wo.)-(a +q^s  *&,„  +a,„)(w,) 


1 42  (  ik  ¥  NO  ] 

2.86  ik^+0j[N0^  +  N0j  ’ 


1  - 1 427  ^  (  k  Y  NO  \ 

2.86  )  [k^+Oj[NO^-i-NOj  ' 


w^ere; 


~  k„„ 


^5(S)/X„  NOj 

k„  +  (‘S's(j)/X„)  1  k^  +  C?j  NO,  +  NO. 


Stage  6 


= o=(a  +e^.  +a„, +a,»)(A'a)-(a +G»„  +a,„ +a„,)(jva) 


n,  A  k„*noi. 


JiWiY (Jf,.XAf02.)Y  o.  \ 

A  n,  A  i:,*N02,  A*-+oJ 
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L42 

2.86 


\ 


NO, 

NO,+NO, 


fi-i-42y. 
V  2.86 


W^iK,){P,)X, 


(  1 

(  NO,  ^ 

j 

[no,+no,) 

>^ere: 


~  ^MO 


^S(6)l-^H 


V 


\r 


K  +  {S,JX„)  ji,  +  O,  A  NO,  +  NO,  J 


NO. 


i^a.) 


Stage  7 

= 0= (a +s,u»  +a<., +a<.,)(w.)-(a +a«  +a<„ +a,„)(  w,) 


Pm(NB) 

( (jr,.)(jvo2,)'] 

f  1 

K  ^NB  J 

1,  K,^N02,  j 

U^+oj 

n-i.42(y.) 
,  2.86 


][(j^«)(^w)(>^,)] 


1.42 

2.86 


*-  1 

r  ivo,  "| 

Uo,  +  AroJ 

^2 


0-i-42y„ 

I  2.86 


^  r 

Fa.(K,)(P,)X„ 


'oh 


yv 


N0,-¥N0J 


V, 
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where: 


k„HSsJX„)^ 


*0/, 


NO, 


k^-^0,){NO,^NO, 


(/^) 


Stages 

o=(a +2«„  +a,..)(Afo,)-(a +a«»  +a,., +a,j,x^o.) 

at 


y.B  A  X,+N02, 


Ak„+0, 


1-1.42(1;) 


2.86 


^][(J^*)(VKf'.)] 


f  Y  w, 

k^+0,]\  NO.  +  NO, 


1-1 42/  ^ 

UK,){P,)X„ 

2.86  y 


Y 


-O* 


V' 


m 


k^-^0,)yNO,  +  NOj 


K 


where: 


=^- 


Y 


k„HSsJX„) 


Y 


A 


k^-^0^ 


NO, 


JK 


NO.-hNO, 


144 


Stage  9 

y  ^ = 0 = (a + +  a,,, + (& + +  a,,, + a«.,)(A'a ) 


V 


(X,,)(N02,) 

A  K,-^N02,  ){k^+0, 


1.42 

2.86 


iVO, 


k^+0,XN0,^N0j 


1-1 427  ^ 


{ 


\( 


A 


NO. 


N0,  +  N0J 


where: 


^S,U  ^ao 


*^S(9)/^W 


KHSsJx„)) 


V 


'■eh 


NO. 


k.,+Oj{NO,  +  NOj 


Stage  10 

ydNO, 

(it 


lo-ft- 


0 = (a> + Qbas  )i^O,)  -  (G, + 0^.  )(NOJ 


Yss  A  K,+N02 


'10 


10  J\  ^«i  ^10  J 


10 
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1.42 

2.86 


V 


V  Ao  y 


mo 

NO,  +  NO,, 


1-IA2Y, 

2.86 


+0, 


10 


m 


10 


NO.  +  NO, 


10 


10  J 


where: 


R  -k 


V 


k^HSsooJ^H)Ak^-^o„ 


NO, 


10 


NO,-^NO,J 


(fa.) 


Nitrite  (NO^*l  Mass  Balances 

(see  mass  balance  equations  for  ititrite  listed  previously  for  Michaelis-Menton 
Model) 

PCOD  Mass  Balances 

(see  mass  balance  equations  for  PCOD  listed  previously  for  Michaelis-Menton 
Model) 
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Oxygen  Consumption 

Stage  6 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  [(/l)(r„  +  l,42(/:,)(X,)](K,)  +  (3.22)(r^.,.,)(FJ  +  (l.U)(r„_,„,)(K.) 


v^iiere: 


o  \  ( 


w 


a 


A=^\-\.A2{Y„) 


'NHtiax) 


_  f  MaKNS)  i 

I  } 


O, 


o,+k. 


'NOiiax) 


(MmjNB) 

y. 


SB  / 


{X^){N02,) 

.  Ks+N02,  , 


a 


V<^6+*, 


Stage  7 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  [(^)(r,  +  1.42(*:,)(jr,)](K,)  +  (3.22)(r^.,„X>;)+(Ul)(r^,,„)(r,) 
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where; 


^S,V  ^ao 


V 


O, 


A  =  1-142(Y„) 


+  k,(P,)(xJ 

y 


o. 


k„  +O7 


’NHtiax) 


’NO^iax) 


V 


f^*l(NS) 

V  ^NB  J\ 


K,+N, 


V 


A 


O. 


o,+K 


\  ^NB  JK 


(X,,)(N02,) 

K^+N02, 


A 


O, 


o,+K 


Stages 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  [(-4)(r.  +  1.42(/rj(jr,)](K.)  +  (3.22)(r^.,.,X)')  +  (l.n)(r^,„,)(K,) 


where: 


^S,V  k^g 


V 


O, 


O. 


^-+0. 


^  =  i-i.42(y„) 
/ 


'Wj(«r) 


^SB  +^* 


Mm(NB) 


V  ^NB  J\ 


{X,,XN02,) 

K,^N02, 


9l. 


A 


0,+kg 


>1 

> 


Stage  9 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  [(^){'^  +  1.42(X,)(^,)](r,)+(3.22)(r,,.,,,)(f;)  +  (l.U)(r„_,„,)(K,) 


where: 


R  =k 

^S.V  "'00 


^5(9)!^ H 


V 


O, 


+  k,(P,)(X,) 


( 


O, 


K^o, 


PmUSSt 


V.  ^NB  J\ 


r 

9  J 


\{ 


O, 


'NOjiax) 


I  A 


v^9  +^«.y 

f 


O, 


^9 


Stage  10 

Dissolved  Oxygen  Consumption  Rate  (Kg  O2  /  hr)  per  stage 

=  [(^X/-,.  +  1.42(/:,)(2f„)](K,.) +(3.22Xr^.,„Xr,.)  +  (1.1  lX'-».,.,Xf:.) 


where: 


R  =k 

^s.v  "mo 


V 


KKHSo,jx„))\k„-^o,J 


a 


10 


^ +  k,{P,,){X,) 


a 


10 


10/ 


y4  =  l-L42(y,  ) 
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V 


K^+N 


\r 


o. 


10 


10  A  *^10 


0,„  +  k 


^m(NB) 


\  ^NB  A 


(X,,)(N02,,) 
K,+N02 


O, 


10 


10  / 


^0,5  +k^ 


